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Abstract 

We fit the 1 Sq A A interaction in the nuclear medium to the masses of the ex- 
perimentally known double— A hypernuclei: ^He, j^Be an< ^ AA-^- We derive this 
effective interaction from OBE Jiilich AA-type potentials and using both Hartree- 
Fock and variational approaches. We find that the inclusion of AA correlations in 
the variational scheme leads to significant differences and a better understanding of 
the dynamical features of the system. We investigate the sensitivity of the binding 
energies and the mesonic decay widths of the above double— A hypernuclei to the 
wAA coupling and the form factor at the oAA vertex. We also use this effective 
interaction to predict binding energies and pionic decay widths of heavier double— A 
hypernuclei, not discovered yet. Finally, we discard the existence of 1 Sq AA bound 
states provided the AA — EN coupling can be neglected. 



PACS: 21.80. +a,21.30.-x,21.10.Dr,21.45+v 

Keywords: single and double— A hypernuclei, AA interaction, Bonn potential, Jiilich 
potential, effective interactions, two-body correlations, mesonic decay. 



Hypernucleus B AA (MeV) AB AA (MeV) 
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10.9 ±0.8 
17.7 ±0.4 
27.5 ±0.7 



4.7 ± 1.0 
4.3 ±0.4 

4.8 ±0.7 



Table 1: Experimental double— A hypernuclei binding energies, -Baa, and values for the quantity 
AB AA defined in Eq. (2). As it is discussed with great detail in p5| ], we have considered that 
the event reported by [7] corresponds to the double— A hypernucleus A \B. We use the notation 

AA Zj - 



1 Introduction 



In the past years a considerable amount of work has been done both in the experimen- 
tal ((l[]-[[ll|]) and the theoretical (||12|-||41|) aspects of the physics of single and double— A 



hypernuclei (for a general overview see the proceedings of the most recent International 
Conferences on Hypernuclear and Strange Particle Physics El). 



Up to now, the experimental community has reported the existence and has measured 
the 1 S'o ground state (angular momentum and spin of the two interacting A's coupled to 
zero, L = S = 0) binding energy of three double— A hypernuclei^]: AA He, \ A Be and AA B 
which energies are reported in Table 0]. Let us consider the hypernucleus AA 2 Z, composed 
of a nuclear core A Z and two bound A particles. The quantity B AA is defined as the total 
binding energy (> 0) of the double— A hypernucleus, and thus it is given by 



B 



AA 



M (ifz) - M ( A Z) - 2m A 



(1) 



where M(- • •) denotes the mass of the system which appears inside the brackets and m A 
(1115.6 MeV) is the A mass. To learn about the nature of the AA interaction, it is usual 
to define the magnitude AB AA as 



AB 



AA 



B 



AA 



2B A , 



where B A (> 0), given by 



M (i +1 z) - M ( A Z) 



m A 



(2) 



(3) 



is the binding energy of a hyperon A in the hypernucleus A +1 Z. Neglecting saturation 
effects, AB AA is suppressed with respect to B AA by one power of the nuclear mass number. 

Data on A-proton scattering constitute an indirect source of information about the 
AA interaction. Although the available hyperon-nucleon scattering data are scarce, this 

1 There has been some confusion about the exact nature of the event reported by Aoki. et. al. [Q. This 
event was identified as or \\B, resulting, respectively, in a repulsive or attractive AA interaction. 

However, the theoretical analysis of Yamamoto et. al. and Dover et. al. |^5| helped to discard the 
possibility of \\Be, and thus it is commonly accepted an attractive nature for the AA interaction. 
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information has been successfully used by the Nijmegen |44j], Jiilich |45|, |46| and 
Tubingen 0, Q groups to determine realistic hyperon-nucleon and thus also some pieces 
of the hyperon-hyperon interactions. The Tubingen model is mainly a constituent quark 
model supplemented in the long and intermediate range part by tt and a exchange; the 
latter is treated as an SU(3) singlet, with a mass of 520 MeV. On the other hand, the 
Nijmegen and Jiilich models are based on meson exchange. In both models most of 
the required intermediate range attraction is provided by the exchange of mesons in the 
scalar-isoscalar channel. 

The Jiilich YN interaction models have been constructed following the same ideas as 
those used in the Bonn NN potential |49|. Thus, these models account for the interme- 
diate range attraction by the exchange of a fictitious scalar-isoscalar meson, a, with a 
mass of about 600 MeV. The a is not treated as a physical particle to which the SU(3) 
relations should be applied, but merely as an effective description of correlated 2n and 
iT^-exchange processes. However, the Nijmegen group views the scalar-isoscalar inter- 
action as generated by genuine scalar meson exchanges (5**, 5, e and k). Then, the 
SU(3) symmetry is used to relate the couplings of the above mesons to the nucleons and 
hyperons. 

Jiilich and Nijmegen models, or some inspired in them, have been extensively used to 
analyze single- A hyper nuclei [ITR O, GIL |3T], p3l 



m. 



made 



issue PL IS, 24, 27 



Because the direct measurement of the AA scattering process is impractical due to 
the lack of targets, the data on AA hypernuclei provide a unique method to learn details 
on the AA interaction in the vacuum. Thus, since the first double— A hypernuclei were 
discovered, suggestions to extract the vacuum AA interaction from these systems were 
Sy]. Since then, many other authors have tried to shed light on this interesting 
, |35], [37j. Near threshold (2m\) , the S (strangeness) = —2 baryon- 
baryon interaction might be described in terms of only two couple channels |]24| : A A and 
EN. For two A hyperons bound in a nuclear medium and because of Pauli-blo eking, 
it is reasonable to think that the ratio of strengths of the AA — > EN — > AA and the 
diagonal AA — > AA (with no EN intermediate states) transitions is suppressed respect 
to the free space case. This is explicitly shown for ®jjie in Ref. |j37| . Thus, the data 
of double— A hypernuclei would mainly probe the free space diagonal AA element of the 
AA — EN potential. In principle, some information about this piece can be extracted from 
the hyperon-nucleon interaction models developed by the Nijmegen, Jiilich and Tubingen 
groups. 



In references [16 



[18| a variational approach and a a-cluster decomposition of the 
considered nuclear cores ( 4 He and 8 Be) are used. However, neither the Nijmegen nor the 



Jiilich or Tubingen models for the AA interaction are used in these references. In Refs. [24 



27], [28|, [37| different Nijmegen AA interactions are considered and G— matrix calculations 
show that the double— A hypernuclei data favor the Nijmegen model D. Despite of the 
success of Jiilich type AN interactions to describe the structure of single-A hypernuclei 
physics [[H], |33|, [SfJ, double— A hypernuclei studies using Jiilich AA one boson exchange 
(OBE) potentials have not been performed yet. 
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Other problems of great interest connected with the AA hypernuclear physics are the 
possible existence of the S = — 2 six quark H dibaryon |5l|, |52fl, and the breaking of SU (3) 
symmetry in baryon-baryon interactions |k|. Six quarks bound states would explore new 
dynamical features of the color interactions different to those accessible by means of the 
study of normal mesons (qq) and baryons {qqq). The existence of AA bound states may 
obstruct the experimental detection of the dibaryon H . On the other hand, some authors 
have shown that the existence of double— A hypernuclei restricts the feasibility of a long- 



living H dibaryon |54 



The experimental community has also invested a lot of effort in the subject, partially 
aimed to discover the H particle (see Refs. |§-|ni). Both at KEK and at BNL facilities, 
cascade particles are being produced by means of the (K~,K + ) reaction and then 
stopped in matter to produce S = — 2 nuclei. At KEK it is planned || to increase 
the number of stopped S~ events by roughly a factor 10 with respect to the experience 
of Refs. 0, U, where the double— A hypernuclei was first discovered. At BNL the 
reaction 

~- A Z - UZ-1) n, (4) 



will be used to produce light and medium double— A hypernuclei 1C, 11 



This work is the first of a series of two where we aim to describe the AA interaction, 
both inside of the nuclear medium and in the vacuum. In this first work we find the effec- 
tive AA potential in the medium. Elsewhere we study the nuclear medium modifications 
of the AA interaction and from this study and the effective potential obtained here we 
will try to extract the AA interaction in the free space, including its spin dependence . 



In this paper, we take the Jiilich model for the interaction between the two A hyperons 
and fit its parameters to the binding energies of the three experimentally known double— A 
hypernuclei. This model does not account for EN intermediate states. Thus, the AA diag- 
onal interaction parameters fitted here, might effectively include some contributions from 
this non-elastic intermediate channel, though we expect them to be reasonably small, 
as discussed above. We use three different approaches: perturbative, Hartree-Fock and 
variational, and find that, though the first two are practically equivalent, the inclusion of 
AA correlations in the variational scheme leads to significant differences and a better un- 
derstanding of the dynamical features of the system. We also study the sensitivity of the 
hypernuclear data to the cutoff masses, used in the Jiilich model for the AA interaction, 
and also examine potentials with dynamical breaking of the SU(3) symmetry. Thus, we 
end up with a whole family of potentials describing the ground state binding energy of 
the three known double— A hypernuclei. With the aim of trying to distinguish between 
the different potentials, we calculate the mesonic decay width of boron double— A hyper- 
nuclei. Finally, we predict binding energies and mesonic widths of heavier hypernuclei, 
not detected yet, we discuss the possible existence of AA bound states and devote a few 
words to relate the free space AA interaction to that found in this work. 

The paper is organized as follows. Our approach to the double— A hypernuclei is 
described in Sect. |2], including the details on the AA interaction. In Sect. |3| our model to 
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the mesonic decay of double— A hypernuclei is discussed. Results are presented in Sect. £|, 
which is split in different subsections: "A-core potentials", "AA interaction: Hartree-Fock 
results", "F 'erturbative approach" , "AA interaction: variational results", "Contribution of 
the 0— exchange" , "Mesonic decay and binding energies of double— A hypernuclei" and 
"Nuclear medium and free space AA interactions. " Finally in Sect. [| we present our 
conclusions. In addition, there is an appendix, where we give some of the needed matrix 
elements within the variational scheme. 



2 Model for the Double— A Hypernuclei 

We approximate the double— A hypernuclei by systems composed by two interacting A's 
moving in the mean field potential created by the nuclear cores (Vaa)- Thus, we solve 



=1,2 



*AA(r!,r$) = 0, (5) 



where Ma and \i& are the nuclear core and the A-core reduced masses respectively, 
^aa^Ij^) is the wave-function of the AA pair and the Vi • V2 piece is the Hughes- 



Eckart (HE) term [56]. The spin dependence in the equation above is implicit in the 
operators and wave-function. The A-nuclear core potential, Vaa, is adjusted to repro- 
duce the binding energies of the corresponding single-A hypernuclei and a a-uo meson 
exchange potential is used for the AA interaction in the medium, Vaa- More details will 
be given in the next subsections. 

In this approximation we neglect the dynamical re-ordering effect in the nuclear core 
due to the presence of the second A and assume that both hyperons move in the same 
mean field as one single hyperon does. Both, the AA interaction and this re-ordering 
of the nuclear core, contribute to A£>aa- However, the latter effect is suppressed with 
respect to the former by one power of the nuclear density, which is the natural parameter 
in all many body quantum theory expansions ]57|]. Thus, we expect the nuclear core 
dynamical re-ordering effects to be of the order of AB\a/A, that is to say around 0.5 
MeV for beryllium and boron and around 1 MeV for helium. However, in view of the 
large binding and incompressibility of an a particle, it is reasonable to think that the 
nuclear-core distortion effects in helium will not be as large as the naive estimate given 
above. Thus, we will assume an uncertainty of the order of 0.5 MeV in all hypernuclei 
due to these effects, which is of the order of the experimental errors of -Baa reported in 
Table |l[ This lack of our model will be translated into an increase of the size of the 
systematic error of our determination of the AA potential in the medium. 

To take properly into account, in the boron region, the effects due to dynamical re- 
ordering in the nuclear core is out of the scope of this paper. Indeed, it would likely require 
the use of variational montecarlo techniques to solve the system of A + 2 interacting 
particles and the use of realistic AW, AA" and AA interactions. The only attempt to 



take partially into account these effects can be found in Refs. [16, 18, pOf for the case of 
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A ^Be. There, the 8 Be nuclear core is described in terms of two interacting a particles 
and the corresponding four-body (aaAA) problem is solved and it is found that nuclear 
core polarization effects change .Baa by about 0.5 MeV, in agreement with our previous 
estimate based in qualitative arguments. This model, however, is difficult to use for any 
other nuclear core that beryllium. 

Note also that the HE term, which naturally appears in our framework, contributes 
to A5aa an d accounts for the variation of the core-kinetic energy due to the presence 
of the second A. As we will see, this piece improves the simultaneous description of an 
extremely light hypernucleus, as helium, and not as light ones, as beryllium or boron. 

Taking into account Fermi statistics for the system of two identical A's, the most 
general wave function in the 1 5'o channel is given by 

^A\(fi,r 2 ) = K[f(r 1 ,r 2 ,r 12 ) + f(r 2 ,ri,r 12 )\ x 5= °, (6) 

where K is a normalization constant, r 12 = \f{ — f 2 \ and x represents the spin-singlet. 

It is also interesting to define the density of probability, P(r 12 ), of finding the two A 
particles at relative distance r 12 . This is given by 

P(r 12 ) = J d 3 n J d^ddn -r 2 \ - r 12 ) |$aa(VI, r^)| 2 . (7) 

We apply both the Hartree-Fock (HF) and variational (VAR) approximations to solve 
Eq. ©. 

In the HF approach, the wave function of the two A's is given by 

*E[(rV2) = Mn)M^)x s=0 , (8) 

<Mr) = R A (r) Foo(r), (9) 

where Yi m (f) are the spherical harmonics and -Ra(^) represents the mono-particle radial 
wave function which is obtained using a self-consistent procedure 



In the VAR approach we allow for a wider family of wave-functions $aa(^1,^), in- 
cluding r 12 — correlations. In analogy with calculations on the atomic three-body prob- 



lem |)8[ - |61[ , we use a series of standard Hylleraas type wave functions [^] to expand 



the wave-function of the two A's system. Thus, our ansatz for the wave function is 



*Ii*(n,f 2 ) 



E Cabc[r\r\ + r\r a 2 ) r^exp {-a(n + r 2 )} 



abc 



x 5=0 , (io) 



with a a real parameter, a, b, c non-negative integer numbers and a < b. Note that, even 
in the case c = 0, where no correlations of the type r 12 are included, this variational 
scheme does not reduce to the HF one because our ansatz of Eq. ([K]) can not be always 
gathered as the product of a function which depends only on the radial coordinate r± times 
the same function of the radial coordinate r 2 . We would like to mention that, though 
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the expected value of the HE term is zero in the HF ground state, it is not when VAR 
wave-functions, with explicit dependence on r 12 , are used. 



For practical purposes, the latter series should be truncated and following the findings 
in atomic physics |59| we choose to fix 



a + b + c< N, 



Thus, the integer number N determines the dimension, Nel, of our basisQ The unknown 
parameters, a and C a b c , are to be determined to give the lowest value of the energy. The 
linear parameters C a bc are obtained by solving a generalized eigenvalue problem p3J, while 



to find the non-linear parameter a one needs to use a standard numerical minimization 
algorithm. We have performed calculations up to iV = 10 for which we estimate that 
the uncertainties in the determination of the binding energies, -Baa, are one order of 
magnitude smaller than their experimental errors. Our ansatz for 3>Xa R an< ^ N = 10 
consists of 161 terms. In atomic physics with such a basis one gets precisions of the order 
of 10" 6 , [ 



61 



much better than those obtained here. That is because the ratio between 
the mass of the orbiting particles and the mass of the nuclear core is much more bigger 
here than in the case of two electrons bound systems. The accuracy could be improved 
by allowing for different exponential behaviors for the radial coordinates r\ and r 2 and 
symmetrizing the resulting wave function 



In the VAR approach, each of the A particles is orbiting around the nuclear core not 
only in the 1=0 wave, as in the HF case. Indeed, if the angular momenta of the two A's 
are coupled to give zero total angular momentum (L = l\ + 12 = 0), any A- nuclear core 
orbital angular momentum is permitted. Of course, the only possibility is that both A's 
have the same orbital angular momentum, h = h = I, with respect to the nuclear core. 

Given a wave function for the two A's system $AA( r i, r 2, r n), the probability, Vi, of 
finding each of the two particles with angular momentum I and coupled to L = is given 
by 



r+00 f+00 i— 1 

Vi = 4tt 2 (2/ + 1) / dr x r\ \ dr*r\ / dp P/(//)$AA(n, r 2 , r X2 
Jo Jo J-i 



:i2) 



where Pi is the Legendre Polynomial of order / and \x is the cosine of the angle between 
the vectors r\ and f*2, being r\i = {r\ + r\ — 2rir2^)3. 



In Refs. [0 and |L8 |, the problem of double— A hypernuclei using a variational ap- 
proach was also addressed. However, the configuration space for the trial wave function 
($^ R ), used in these references to describe the AA relative motion, is much more reduced 
than the one used in this work. 



2 A^el is given by 



iV EL = 



i (2iV 3 + 157V 2 + 347V + 24) N even 
^(2iV 3 + 15iV 2 + 34iV + 21) N odd 
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With a basis of the size of the one used here, one needs to compute more than twenty 
five thousand matrix elements of the Hamiltonian, for each value of a. Furthermore, be- 
cause of the non-orthogonality of the chosen basis there exist large numerical cancelations 
in the linear parameters (C a b c ) minimization process and it is necessary to have analytical 
expressions for the matrix elements. These can be found in the Appendix. 



2.1 Model for the AA Interaction 

Given the quantum numbers of the A particle, the lightest carriers of the strong interaction 
between two A's, in the scalar-isoscalar and vector-isoscalar channels, are the a (I = 
0, J p = + ) and uj (I = 0, J p = 1~) mesons, respectively. As it was mentioned in 
the introduction, the a is not treated as a physical particle, but merely as an effective 
description of correlated 2tt and .fT.K'-exchange processes, as it is done in the context of 
the Jiilich models for YN interaction H5L Bfl. The 77 and 7/ meson exchanges lead to a 



AA potentials in the pseudoscalar-isoscalar channel which contributions are negligible at 
low momentum transfers R, and are not considered in the context of the Jiilich models for 



YN interaction |45|. [46| . Thus, we will not consider those in this work. 



In the spirit of the Bonn- Jiilich potentials, we should also study the exchange of the 
<p (I = 0, J p = 1 _ ) meson. In Ref. f35|, where the Nijmegen model D contributions of 
different mesons are studied and compared for the 1 S'o AA system, it is shown that the 
contribution next in importance to that of the a and uj mesons is the <fi— exchange, being 
its contribution attractive. By taking the so-called "ideal" mixing ang le (tan0„ = l/y/2), 



the <fi meson comes out as a pure ss state and thus in Refs. [45], [46| a vanishing cj)NN 
coupling is required. This provides a relation between the singlet and octet couplings, 
which determines the 0AA couplings in terms of the uAA ones. 

The mass of the <fi meson is significantly larger (m^ = 1019.41 MeV) than those of 
the a and uj mesons. Then its contribution will be relevant at very short distances, where 
any potential model based on meson exchanges suffers of some uncertainties which are 
translated into the inclusion of form-factors and uncertainties in the precise values and 
shapes of them. To explore the dynamics at short distances one needs high momentum 
transfers, thus though one might expect the 0— exchange to be relevant in scattering 
processes at large energies and scattering angles, it is also reasonable to think that such 
contribution is much less relevant to study the dynamics of a AA pair below threshold. 
Furthermore, because the meson does not couple to nucleons within the Bonn- Jiilich 
model, its contribution was not taken into account in their phenomenological studies of the 
NN and YN scattering data. Thus, we count with any empirical determination neither 
of the <f)NN nor of the 0AA form-factors. We have checked that the importance of the 
0— exchange in the 1 So A A system at threshold is small and, depends strongly on the used 
form-factor, as one could expect. Taking as reference the cutoff used for the uj— exchange, 
we find that small changes in the cutoff lead to big changes in the small role played by the 
(f) meson in double— A hypernuclei. Given that we count with only three piece of data, that 



This is explicitly shown in Refs. |35L H also for the Nijmegen model D. 



7 



we expect the 0— exchange contribution to be small and furthermore is not completely 
determined by that of the uj— exchange and thus its inclusion would require additional 
free parameters, we will not include the 0— exchange from the beginning. We will try to 
keep our model for the effective AA interaction in the medium reasonably simple, and thus 
we will construct it in terms only of a— and uj— exchanges. In any case, in Subsect. ET5 



we will discuss qualitatively and quantitatively the influence of the 0— exchange in the 
results presented through the paper. 

Apart from single meson exchanges, in principle, one should also include the AA — 
EN coupling, as mentioned in the Introduction. For simplicity, we will first ignore this 
possibility. The possible contribution of this coupled channel and that of the <fi meson 
exchange (and that of heavier mesons, not explicitly mentioned in this discussion) will be 
in any case accounted for in the obtained effective coupling constants of the considered a 
and uj exchanges. We will comment on this when we try to link the effective AA interaction 
obtained in this work, from the double— A hypernuclei data, to the interaction in the free 
space. 

The couplings of the a and uj mesons to the A are described by the following interaction 
Lagrangians |36|, f4"5| : 



AxAA = 9aAA^A{x)(f>a{x)^A{x), ( 13 ) 

CuAA = gu,AA*A(xhMx)**{x) + ^^ A (x)a^ A (x) {d»<%(x) - d v ^(x)) (14) 

4m N 

where we have used uin = 938.926 MeV for the nucleon mass. \&a, 4><t an d 4>uj are the 
A, er and uj fields, 7 M are the Dirac matrices and a^ v = ipy^, 7 I/ ]/2. The parameter /^aa 
accounts for the tensorial part of the u;AA coupling. 

The AA OBE potential in configuration space is obtained via Fourier transform from 
the AA Feynman amplitudes 4 obtained from the above Lagrangians and reads [49 



V a (m a ,r) 



V u (m u ,r) 



+ 



+ 



n 2 

9aAA 

4:71 



1 _ I ( m °\ 2 
4 \mj 



Y(m n r) 



4m\ 

n 2 

Am 



V 2 Y(m a r) + Y(m a r)V 2 



-Zi{m a r)L ■ S 



(15) 



1 + 



1 /m^A 2 

2 \mj 



Y(m,,r) 



4m\ 



V 2 y(m w r) +F(m w r)V 2 



6 \m\J 



12 



2 4tt m N {\m A / 3 \m A / 



m^r at ■ a 2 



4 A non-relativistic expansion of the amplitudes is performed, keeping only up to quadratic terms in 
momenta over baryon masses. 
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- AZ 1 (m w r)L ■ S - -Z(m w r)Si 2 



(fi-a 2 - -2(mj)Si 2 



Y(x) 



e~ x „, 1 d 2 L 
, V 2 = -— r- 



(16) 
(17) 



x r dr 2 r z 

where m a uj\ is the a (u) meson mass for which we take 550 (782.6) MeV. The operators 
a, S, L and Su are the Pauli matrices, total spin and angular momentum and the second 
order rank spin tensor operator respectively. Finally, the functions Z, Z\ can be found 

g. 



in 



The potentials should also contain form factors describing the extended hadron struc- 
ture. Although they in general depend on all the three four-momenta involved at the 
vertex, they are usually parameterized in a simple form depending only on the four- 
momentum, q, of the exchanged meson. Taking a monopolar form and neglecting the 
energy transfer dependence at the vertex [49], [45[ , we use form factors of the type 

F Q (q) = ^ AA ~ m " , (If 



A 2 -I- a 2 ' 



where the a symbol stands for a or uj. 

The use of the given form factor at each vertex leads to the following extended ex- 
pressions for the potentials: 

{A 2 aAA -m 2 a ) dV Q (z,r)\ 



V a (r) = V a (m a} r) - < V a (z,r) 



2A 



aAA 



dz 



(19) 



Z=A a AA 

Furthermore, for the 1 Sr t channel, in which we will be interested in this work, one can 



replace 



01 ■ CT 2 
L -S 

Sl2 



-3. 



0. 



Thus, finally and after including the above substitutions, the AA potential reads 



V, 



L=S=0 ( 
AA \ 



L=S=0 



(r) + Vf= s =°(r) 



where 



L=S=0/ 



n 2 

9aAA 
47T 



m a I Y(cr, r) + 



VY(o\r 



2m\ 



V + Y(<x,r)V z 



(20) 
(21) 
(22) 

(23) 
(24) 



V. 



L=S=0i 



Air 



jAA 



(91aa - sIaa) ( a Iaa - m 



2 \2 



-A„aa r 



m:. 



w 1 ^luAA 



ojAA 



2m\ 



(\7Y(cu,r)) L - V + Y(w,r)V 2 



(25) 
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Vertex 


9a/ 


fa/V^ 


A a (GeV) 


uAA 


2.981 


-2.796 


2 


aAA 


2.138 




1 



Table 2: Coupling constants g a , f a and cutoff masses A a found in Ref. [46] for model A 



The subindex L means that the operator V only acts on the function Y(a,r). On the 
other hand, g^AA and Y(a,r) are given by 

6- 2 n 2 1 (m^ 2 ( 3g* AA m A /™a\ 2 , 2 \ 

2 \mA/ \ 2 \m N J j 

Y(a,r) = Y(m a r) -\l + 7 ^(Al AA - ml)} ^Y(A aAA r). (27) 



The phenomenological Lagrangians of Eqs. (|13D - (|TJ) and form factors of the type 
of Eq. ([Tj|) have been used in Refs. ||| to describe AN elastic scattering dataQ. The 



interactions of Ref. ]1H] are obtained from those in Ref. [|5|| by neglecting their small time 
dependence and performing again a fit to the same set of data. The NN vertices are 
taken from the Bonn potential [39"]. The model A of Ref. leads to a AA potential 



of the same form of that of Eqs. fljjD - (0) and in Table ^| we give the values of the 



coupling constants and cutoff masses found in this reference. Note that in Ref. ||46|| , the 
strength of the cuAA coupling is fixed by means of the SU(6) symmetry (g^AA = \g U NN 
and /ojaa = \fuNN — \IpNn) and thus g a AA, A ct aa and A^aa are the only free parameters 
adjusted to the empirical data. In that reference is also pointed out that the fit has little 
sensitivity to the cutoff mass A ctA a- 

In principle, in this work the free parameters would also be the coupling constants of 
the a (<7ctaa) and to (g^Ak and / w aa ) mesons to the hyperon A and the cutoffs (A ct aa, A^aa) 
of the corresponding form factors. However, due to the limited set of data available to this 
study, we fix A^aa = 2 GeV and / w aa = —2.796, as in Table |2|. To check the sensitivity 
of the hypernuclear data to g u Ah and A ct aa, we consider different families of potentials 
obtained by using different values for the ratio g^Ak/ 9wNN above and below the SU(6) 
prediction, 2/3, and taking for each value of the former ratio, two different values (1 and 2 
GeV) of the cutoff A ct aa- Thus, for each potential we are left with just one free parameter, 
<?<tAA, which is obtained from a best x 2 fit to the values of -Baa reported in Table [I]. 

As we will see, there is a strong correlation between the parameters g a AA and AoAa, 
and the quantity g a AA X (1 — m 2 ,/ ' ^oAa) remains constant within approximately 10%. This 
is due to the fact that for bound A particles in the double— A hypernucleus, the average 
momentum transfered in the AA vertex is much smaller than the typical values (~ 1 

5 In these references, it is not only studied the AiV — > AiV process, but all hyperon-nucleon reaction 
channels measured experimentally (AiV —> AiV, EiV — > AiV' and SiV — > S'iV') using a coupled channel 
formalism. 



10 



GeV) of the cutoff mass. In these circumstances the form factor at the aAA vertex can 
be approximated by 



2 / -»2 

m a I q 



A 2 \ A 

A oAA \ A cr AA, 



Similar conclusions can be drawn in the case of g^AA an d A^aa- 

Finally, we would like to mention that the AA interaction determined in this work 
corresponds to an effective interaction in the nuclear medium and it is not directly com- 
parable with that of Ref. [iB|. We will come back to this point in Subsect. W7J\ 



2.2 Model for the A-Core Interaction 

The detailed study of the A-nuclear core dynamics is not the aim of this paper. However, 
we need to check the sensitivity of the parameters of the AA interaction, obtained from 
an overall fit to the available -Baa data, to the used A-nuclear core potential. Thus, we 
have considered different A-nucleus potentials suggested in the literature, which have been 
adjusted to reproduce the binding energies Ba of the single-A hypernuclei ground states. 

In this subsection we will give details of the different A-core potentials used through 
this work. All except for one are obtained from models to describe the dynamics of the AA" 
system. To obtain a A-nucleus potential (Vaa) from a AA" interaction (Van) we fold the 
latter with the corresponding nuclear density of the core. For closed-shell nuclear cores 
and A particles in the ground state (Is a), all the couplings involving the total angular 
momentum L or spin couplings of the type &[ ■ a 2 or S12 in the Van potential do not 
contribute, and thus for this particular scenario the Vaa potential is given by 

Vaa(t) = I d 3 r' Pc (\f- r ' \)V AN (r'), (29) 



where p c corresponds to the density for the centers of nucleons which is obtained from the 
charge density by taking into account the finite size of the nucleon (for details see Sect. 4 
of Ref. H). 



A considerable amount of work has been devoted to the study of the mesonic decay of 
the hypernucleus A He as a mean to investigate the repulsive part of the AA" interaction. 
Then, different types of A- 4 He potentials have been suggested in the literature |22|, p9|| . 
Among them we have selected two potentials, with (YNG) and without (ORG) a repulsive 
part in the elementary AA" interaction. These potentials are determined by the following 
AA" interactions: 

Van(t) ORG = -38.19 e-^ L034 ) 2 MeV, (30) 

V AN (r) YNG = [919 e"^ - 5 ) 2 - 206.54 e~^ 2 - 9.62 e -('A- B > a ] MeV. (31) 

In all cases r is given in fermis. The parameters of the above interactions have been 
adjusted to reproduce the ground state binding energy of the A He f22| . 

To perform a simultaneous analysis of the three double— A hypernuclei we are inter- 
ested in, we have also considered other two A-core potentials. 
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In the spirit of the previous subsection for the AA dynamics, the first A-core potential 
(SW) considered here, is obtained from an effective a—uj exchange model to describe 
the KN interaction given by 

V AN( m <r,r) = - 9aAA9aNN m^Ylm^r), (32) 

47T 

V\N( m u, r) = — < g^wgwNN - -. 

47r [ 4 \m A m N J 

X ( guAAOuNN —(juAAfuNN ~ <7wi\wA>AA ) \ Y(m w r). (33) 

V m N J J 

The potential above is a simplified version of the more general one described in 
Eqs. (38) - (39) of Ref. Firstly, as we will mention below, we use monopolar 



type form factors instead of Gaussian type form factors as it is done there. Secondly, 



we have neglected all terms of order O (Jm mcson /m baryon ) J and also all terms of 

order O (^meson( m A — m Af)/ m baryon)- Thirdly, we have neglected all type of spin 
terms which will not contribute for closed-shell nuclear cores and A particles in 
the ground state {Is 2). Finally, we have neglected non-local terms of the potential 
which would be proportional to the A momentum squared over the product of the 
A and nucleon masses (• ■ • (V 2 0A(^*))/4mA^Ar). For a bound particle, one expects 
these contributions to be quite small. 

We do not consider potentials due to the exchange of strange mesons (K, K*, • • •), 
which though relevant to describe the AA^ scattering process [f43|, |4"6|| , contribute to 
the binding energy only as Fock terms and thus are suppressed at least by a 1/A 
factor. In Ref. ]TJ| these Fock contributions were evaluated and found to be small. 

To take into account the finite size of the baryons involved at the vertices, monopolar 
form factors are used here as well. In coordinate space the inclusion of form factors 



is implemented by means of the substitution [4G 



A 2 2 

VMr) = VZ N (m a ,r)- « NN ~ ™ a ^v(A aAA ,r) 



A 2 _ A 2 

Iv aNN Iv aAA 



i" AA ?2 a V£ N (A aNN ,r), (34) 

A oJVV A oAA 



where the a symbol stands for o or u. Thus, finally the AA^ potential reads 

V^(r) = VUr) + VUr)- (35) 



In principle, the couplings g w ^A, fi'oAA an d /Laa introduced in Eqs. (JT3) - (|16|) 
should coincide with those appearing in Eqs. Q3"2"D - (|3"BD, and denoted there with 
an overline. However, these couplings might be significantly different, because in 
both cases we are dealing with effective interactions (AA or AA^ ) which parameters 
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Vertex 


ga/V^TT 


/o/v^F Aa (GeV) 


uNN 


4.472 


1.5 


aNN 


2.385 


1.7 



Table 3: Coupling constants g a , / a and cutoff masses A a for the AW vertices. Values are taken 
from Ref. II. 



also account for some contributions (0— exchange, intermediate states: 'EN, UN, 
etc..) not explicitly included. Furthermore, these effective interactions are affected 
by renormalization effects due to the nuclear medium, which could be significantly 
different for each interaction. On top of that, as we mentioned above, some non- 
localities and Fock contributions have been neglected in the treatment of the A-core 
interaction, and their effect, though presumably small, could affect the value of 
the effective couplings in Eqs. (|32] ) - (|33|) . Thus, to avoid confusions with the 
AA couplings, we denote with an overline those couplings which appear in the AN 
potential. 

In Eqs. (RBI) - (51) we fix the uAA coupling parameters to the values given above 



in Table 0. The uNN and aNN coupling parameters are taken from the Bonn 
potential [[l9j and compiled here in Table |3|. Finally, for each hypernucleus (^ +1 Z) 
we fit g a \A to reproduce the ground state binding energy, using two values, 1 (SW1) 
and 2 GeV (SW2), for the cutoff parameter A^aa, as we discussed in the previous 
subsection. 

The second A-nuclear core potential (BOY) considered, is phenomenological and it 
is not based on any model for the AN interaction. It was suggested long time ago 
for medium nuclei by A. Bouyssy [|13j. It has only one parameter, Vq, which we 
adjust for each hypernucleus (f +1 Z) to reproduce the ground state binding energy. 
The potential reads: 

v boy ( s _ v o 

^AA V > l _|_ e (r-R)/a ' 

R = LlAs fin, a = 0.6fm. (36) 

A similar potential was also suggested in Ref. |Tj|. There, it is also shown that, 
despite of not including a spin— orbit part, these type of potentials give also rea- 
sonable descriptions of non s-wave A bound states, in agreement with the previous 
findings of Ref. Il2| . 



3 Mesonic Decay of Double— A Hypernuclei 



The mesonic decay has been computed following the method exposed in Refs. PO] and PD 



which uses shell-model baryon wave functions and distorted pion waves. For light nuclei 
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in the p shell, as boron, the low lying state structure might not be correctly represented 
simply by unoccupied single-particle orbitals, and a residual interaction of the Cohen- 
Kurath type |J7| has been added. Reviews on the subject can be found in Refs. |33| and 



41] . We compute the decay widths corresponding to the following processes 



AA n 

A+2 7 
AA L 

AA L 

A+2 7 
AA L 



A+2r 



7T 



+ n + 7r 



o 



9* 



(37) 
(38) 
(39) 
(40) 



where b denotes that the remaining A is in the ground state and the outgoing nucleon is 
in an unoccupied bound state of the daughter hypernucleus. On the other hand gs means 
that the daughter hypernucleus is left on its ground state. In the two last reactions the 
outgoing nucleon goes to the continuum. 

We use a model in which the pionic decay is produced by a one-body operator 



SH 



AnN 



-Gml 



P 

s-{—) 



a ■ q 



T 



(41) 



where (Gm 2 ) 2 /8n = 1.945 10 -15 , the constants 5* and P are equal to 1.06 and 0.527 
respectively and is the pion mass (139.6 or 135.0 MeV for 7r~ or 7T°). Finally, q is 
the momentum of the outgoing pion and the Pauli matrices a and r A act on the spin and 
isospin Hilbert spaces respectively. The r A operator in Eq. (f4"T|) implements the AT = 1/2 
rule by means of which the rate of A — > ir~p is twice as large as that of A — > 7T°n. 

The vacuum A decay width is readily evaluated and leads for proton or neutron decay 

to 



4«) 

free 



(Gml) 



m N q C m 




(J (a) 

4ir ttia 
\ 1/2 (m A ,m 2 N: ml) 
2mA 

x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, 



(42) 



(43) 
(44) 

the pion momentum in the 



X(x,y,z) 

with the isospin coefficients =4 and = 2, and q, 
center of mass frame. One can see from Eq. (|42"D that the parity violating term, S, is the 
dominant one in the decay. 

To illustrate the main ingredients entering in the decay, we reproduce below the decay 
width for any of the processes of Eqs. ( |3"7| ) - ([H]), i n the simple shell-model case where 
the spin-orbit splitting of the nuclear-core levels is not taken into account: 

d 3 q 1 



Y ( Q ) 



C (a) £ J 



X 



+ 



S 2 



J d 3 x 1 



(2tt) 3 2u(q) 



2vr(5(m A - B AA + B A - u(q) - E N )(Gm 2 r 



2\2 



<fi ) {q,xi)*^* N (x l ] 



' p \ 2 



d 3 xi 



d 3 X 2 $AA{Xl, X 2 ) <f* A (X 2 ) 



Vtifl ) {q,x 1 )* V * N (x 1 ] 



(45) 
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where tp^ and En are the wave function and energy of the outgoing nucleon in the A 
decay, $aa is the wave-function of the AA pair, ip^ is the A wave function in the daughter 
hypernucleus, u>(q) the pion energy, and the sum over N runs over the unoccupied nuclear 
orbitals given by n and / since spin sums are already performed. In Eq. (|45|) the sums over 
N are over proton or neutron orbitals according to a. Note that if the AA interaction and 
the HE term were neglected and if B\\ were replaced by 2B^ in the energy conservation 
delta, then Eq. ( pE5| ) would yield a width of twice the corresponding one for the decay of 
a A in a single-A-hypernucleus. 

The corresponding expression for the decay width when the spin-orbit splitting of the 
nuclear-core levels is considered, can be easily deduced from Eqs. (6)-(15) of Ref. [[20| . 



The pion wave function (cp(. \q,x)*) as a block corresponds to an incoming solution 
of the Klein Gordon equation, 



-V + ml + 2uj(q)V opt { 



x 



^-\q,xy = (u(q) - V c (x)) 2 ^-\q,x)\ (46) 



with Vc(x) the Coulomb potential created by the nucleus considering finite size and 
vacuum polarization effects and V opt (x) the optical potential which describes the 7r-nucleus 
interaction. This potential has been developed microscopically and it is presented in detail 



in Refs. |68], |>9| . It contains the ordinary lowest order optical potential pieces constructed 
from the s- and p-wave irN amplitudes. In addition second order terms in both s- and 
p-waves, responsible for pion absorption, are also considered. Standard corrections, as 
second-order Pauli re-scattering term, ATT term, Lorentz-Lorenz effect and long and 
short range nuclear correlations, are also taken into account. This theoretical potential 
reproduces fairly well the data of pionic atoms (binding energies and strong absorption 
widths) p8l and low energy 7r-nucleus scattering [p9(| . 



The A wave function in the daughter hypernucleus, y?A, is computed using the A- 
core potentials described in Subsect. ^2| . For the nucleons we have used the following 
potential J7D| 

-5° MeV 

^> = 1+exp[(r - fl)/a] - < 47 > 

with R = 1.25A 1 / 3 fm, a = 0.65 fm, which provides a fair reproduction of the nuclear 
levels for the average energy of major shells, as well as realistic nucleon wave functions. 

Because, the strength of the spin-orbit force is small in the boron region (about 4 
MeV), we will neglect it, and to calculate the mesonic decay of the hypernucleus 
we will use the potential of Eq. (f47|) and Cohen-Kurath spectroscopic factors [|7l[] for 
the lp— shell. Those are obtained from effective interactions for the lp— shell deduced in 
Ref. |[7j by fitting different nuclear energy levels. Thus, the contribution to the mesonic 
width of the processes where the outgoing nucleon is trapped in the lp— shell is computed 
as follows. The initial state consists of seven nucleons in the lp— configuration^] plus two 

6 The lsi — shell is full, and it behaves, in a very good approximation, as an inert core 
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A hyperons coupled to L = = S, whereas the final state consists of eight nucleons in the 
lp— configuration plus a A hyperon and a pion. The transition operator annihilates a A 
hyperon and creates a nucleon in the lp— configuration. To compute the nuclear matrix 
elements involved in the mesonic decay, one needs the angular momentum-isospin reduced 
matrix element of the nucleon creation operator, aL between fully antisymmetric states 
of n and n — 1 nucleons in the angular momentum-isospin configuration p and coupled 
respectively to T R and Tj angular momentum and isospin quantum numbers. This is 
given by [731, 



(p n r R \\\al\\\p n - l ri) 

^ny/2T R + 1 



where the reduced matrix element is defined by Eq. (A.3.al7) of Ref. [73| and (...|}...) are 



the coefficients of fractional parentage. They are related to the spectroscopic factors, S, 
by 

S{T R ;Ti,p) = n(p"r i? |}p"- 1 r,) 2 . (49) 

After a little of Racah-algebra, one finds that the averaged modulus squared of the 
nuclear matrix element (Ai) which determines the lp— shell contribution to the decay 
processes of Eqs. (j37l) - fl38|) is given by 



i f 



E ( 2 j. T^7^~ i) C ( Ti > \ T R\ T ^ t R ~~ r i' t r) 2 S{Jri t Ri "a; Ji, T ii gs, j] 



I 1 ~ 1 \ 2 

22 (A-mA s \8H A7TN \Njm,-(T R -Ti) ®nq(a R )) 



X 

m,m As 



(50) 



where Jj = 3/2, Tj = 1/2 and Tj = —1/2 are the total angular momentum, isospin and 
third isospin component quantum numbers of the initial nuclear core ( n B), J R ,T R and 
t r = for tt~ decay or —1 for tt° decay, are the corresponding quantum numbers for 
the final nuclear core. In addition gs and a R stand for the ground state energy of the 
n B core and the excitation energy of the final cores 12 B* or 12 C* and C is a Clebsch- 
Gordan coefficient. On the other hand j, m and t r — n determines the wave function of 
the outgoing nucleon after the decay. The sum over the "R" quantum numbers gives us 
the sum over p— shell excited states in the final nuclear core. For 7r~ (tt°) decay, we sum 
the contribution of the eighteen (eight) excited states which spectroscopic factors and 
energies are given m Table lQ (label: "Stripping for target A = 11 (§, |)") of Ref. 



Finally, \irq(a R )) and |AmA s ) stand for the pion wave function in the continuum, being 
its momentum determined by energy conservation, and the spectator A wave function, 
being its spatial wave-function determined by the projection 

7 Note that in this table the origin of energies is that of the 12 C ground state. 
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(^a(x 2 )\^aa(xi,x 2 )) = / d 3 x 2 ® AA (x 1 , x 2 )ip A (x 2 ). 



(51) 



The matrix element of the interaction 5H AnA f in Eq. (pOj) is of the same type of that 
given in Eqs. (6)- (15) of Ref. ED . 



Further details of the calculation of the mesonic decay width as the treatment of the 
outgoing nucleons in the continuum, the correct energy balance in the reaction, the correct 
treatment of the quasielastic collisions of the outgoing pion, the procedure to perform the 
d 3 q and dfl\ integrations, etc., can be found in Ref. |30J . 



Experimentally, what can be observed are the inclusive processes 

a+2 Z ^ x + A + n , (52) 
IfZ -> X + A + 7T-. (53) 

The main contribution to these processes is given by the exclusive ones shown in Eqs. (|3"7|) - 
( flUD . This can be seen by looking at the overlap integral defined in Eq. ([51]) , which appears 
in Eq. (fl5}). If one uses HF wave functions for the AA pair and the pure spectator approx- 
imation, in which the monoparticle wave function used to construct 3>aa coincides with 
that of the A ground state in the daughter hypernucleus, then, due to the orthogonality 
of the A wave functions, only those processes where the non-decaying A in the daughter 
hypernucleus remains in the ground state will contribute. Even if these approximations 
are not used, we have checked that those processes included in Eqs. (ET) - (Ek3) and not 



included in Eqs. ( j3?|) - (f40|) represent less than 2-3% of the total. Therefore, in what 
follows we will approximate the inclusive double— A hypernucleus mesonic decay widths 
(to 7T° or to 7r~) by those corresponding to the processes specified in Eqs. (p7|) - (f40|). 



4 Results 



Before we start discussing the results concerning double— A hypernuclei, we show briefly 
the results for the A-core potentials described in the Subsect. |2~2| . 



4.1 A-Core Potentials 



We have studied not only the three experimentally observed double— A hypernuclei. To 
study the A dependence of AB AA and the mesonic decay width, we have also consid- 
ered the nuclear core closed-shell hypernuclei A \C<i, ^Zr and ^Pb. Unfortunately, the 
binding energies of the corresponding single-A hypernuclei ( 4 \Ca, 9 \Zr and 209 Pb ) are 
not known and we have had to approximate them by the binding energies of the closest 
hypernuclei known experimentally. Thus, we have used the experimental binding energies 
corresponding to 4 A Ca, S9 A Y and 20 |Pb. In Table [| the binding energies, A-core reduced 
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Hypernucleus B\ xv [MeV] density A R (c) [fm] g(g[fm]) [i A (MeV) 





3.12±0.02 0, 


HO 


4 


1.358 





858.6 


a Be 


6.71±0.04 


HO 


8 


1.77 


0.631 


970.4 


IB 


11.37±0.06 [0 


HO 


11 


1.69 


0.811 


1006.1 


\ Ca 


18. 7± 1.1 § 












a CcL 




2pF 


40 


3.51 


0.563 


1083.17 


89 v 
A 1 


22.0± 0.5 § 












«i Zr 




2pF 


90 


4.84 


0.55 


1100.95 


20 A Pb 


26. 5± 0.5 § 












20 A Pb 




2pF 


208 


6.624 


0.549 


1109.21 



Table 4: Binding energies, charge densities (JF^]), nuclear mass numbers, and A-core reduced 
masses for several single-A hypernuclei used through this work. 2pF stands for a two parameter 
Fermi density, p2pF{f) = Po/ fl + e^ r ~ c ^ a ^j, and HO stands for an harmonic oscillator density, 

Pno( r ) = po(l + a (r/R) 2 ^j e~( r l R } 2 . The density of nucleon centers, which appears in Eq. (pS|), 
is obtained from the charge density by taking into account the finite size of the nucleon (for 
details see Sect. 4 of Ref. |]6q1 ). 

masses and details of the charge densities for all single-A hypernuclei considered in this 
work can be found. 

In Table [5], we show the binding energies and mean squared radius obtained for the 
ground state of different single-A hypernuclei by using the A-core potentials described in 



Subsect. pT2|. In this table we also give the parameters V and g a \A fitted to the ground 



state binding energy for each hypernucleus. We would like to make a few remarks: 



YNG and ORG potentials were only adjusted in Ref. |22| for A He. Thus, we will 
restrict their use only to the case of helium. 



In the original work of Ref. |T3[ an overall fit to light-medium hypernuclei was 
performed which provided a value of Vo = —32 ±2 MeV, which is in good agreement 
with the values shown in Table [5[ 

One might think that for the case of the o — uj interaction with a cutoff A ct aa = 
1 GeV, the fitted parameter g a AA/V^ should coincide with the value found in 



Ref. |[46|| , 2.138, for hyperon-nucleon scattering in the vacuum. We find significantly 
different values, 3.1-3.5, which may be due to the fact that we are dealing with 
a medium interaction, that we have neglected Fock terms and some non-localities 
which may be significant for light and heavy nuclei respectively, that we are using 
a A-core potential designed for closed-shell cores, for some hypernuclei where this 
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is not the case, and finally that we have not explicitly considered here the channel^ 
AN — > TjN which is included in the analysis of Ref. ||46|| . To trace back the origin 
of this discrepancy is even more complicated when one realizes that, in both this 



work and that of Ref. pig] , for a given set of cu— couplings, what is determined is the 
product of QoKKgoNN rather that the g a AA~ coupling alone. Hence, the observed 
discrepancy could also be due to differences between the aNN coupling in the 
vacuum and inside of the nuclear medium. 



In contrast to the findings of Ref. [0, our fits depend appreciably on the cutoff 
parameter A a ^A- Indeed we find that the quantity g a \A X (1 — m 2 a /K 2 aKh ) remains 
constant within approximately 2%, as we expected from our discussion in Eq. (|28|). 
Despite of the fact that the latter quantity remains almost constant, for very light 
hypernucleus as ^He, the short-distance behavior of the potential depends strongly 
on the specific value of the cutoff A^aa, as can be seen in Fig. (E[ 

Because of the folding in Eq. (^), the repulsive or attractive character at short 
distances of the AA" interaction can be appreciably seen in the A-core potential 
only for very light hypernuclei as ^He. This is clearly shown in Fig. [I], where one 
can also see that for medium and heavy hypernuclei none of the A-core potentials 
provide repulsion at short distances despite of the fact that some of them are being 
constructed out of AA" potentials with repulsive cores. 

For a — uj type potentials, there is a strong cancelation between both components 
of the interaction which is translated into a great sensitivity of the binding energy 
to very small changes of the coupling constants, as the extremely small statistical 
errors of the parameter g a AA in Table ^| indicate. 



4.2 AA Interaction: HF Results 

In Table ||| we show, within the HF scheme, the parameter g a AA obtained from a best fit 
to the available B^a experimental data (Table [I]). The value of the u;AA coupling is kept 
fixed to g^AA = tyuNN /3 and we consider several A-core potentials and cutoff parameters 
AoAA- The x 2 by degree of freedom for each fit is also shown. YNG and ORG potentials 
have been considered only for helium. 

The fits using the ORG for helium and BOY for beryllium and boron A— core in- 
teractions, provide the best x 2 /d°f (~ 0.1) among all A— core interactions considered. 
Anyhow, most of the fits presented in Table |6] with different A— core interactions are sta- 
tistically acceptable, and in all the cases the differences in g a AA from the different fits 
are of the order of the statistical errors. Indeed, the fitted values of g a AA only depends 
significantly on the value of the cutoff parameter A ct aa- Thus, the ambiguities in the 

8 Though relevant far beyond the AN threshold, we expect the contribution of this new channel to 
be small for A-bound states because of Pauli-blocking |}2| . In any case it will contribute to the above- 
mentioned difference between the <? ct aa couplings. 
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YNG ORG BOY SW1 SW2 

A.aa [GeV] 1 2 

Hypernucleus 



iHe 


,B A [MeV] 3.12 

/ 2\l/2r f i o 16 

\ * / J- -Li- J. O • -L \J 

par: -V or ^ A A 


3.10 3.12(2) 
2 77 3 02(1) 

— .11 UiUil 1 1 

29.95(6) 


3.12(2) 
3.12(1) 
3.4916(11) 


3.12(2) 
3.00(1) 
2.5662(7) 




fi A [MeV] 
/ 2\l/2r f i 

par: -V or g aAA 


6.71(4) 
2.46(1) 
28.02(7) 


6.71(4) 
2.61(1) 
3.4389(16) 


6.71(4) 
2 533(5) 
2.5436(12) 


IB 


B A [MeV] 

( r 2^1/2[ fm ] 

par: -Vq or ^ A A 


11.37(6) 
2.216(3) 
31.98(10) 


11.37(6) 
2.259(4) 
3.3938 (14) 


11.37(6) 
2.183(3) 
2.506(1) 


4 A Ca 


5 A [MeV] 

( r 2)l/2[ fm ] 

par: -V or g aAA 


18.7 ± 1.1 
2.47(3) 
29.7 ± 1.3 


18.7 ± 1.1 

2.21(3) 

3.194(16) 


18.7 ± 1.1 

2.26(3) 

2.390(12) 


9 A Zr 


B A [MeV] 

( r 2)l/2[ fm ] 

par: -Vq or g CTA A 


22.0(5) 
2.91(1) 
28.9(5) 


22.0(5) 

2.562(7) 

3.129(7) 


22.0(5) 

2.645(8) 

2.353(5) 


20 A Pb 


fi A [MeV] 

( r 2)l/2[ fm ] 

par: -Vq or g aAA 


26.5(5) 

3.588(7) 

30.7(5) 


26.5(5) 

3.341(5) 

3.145(7) 


26.5(5) 

3.428(5) 

2.371(5) 



Table 5: Binding energies and A-mean squared radius obtained for the ground state of different 
single- A hypernuclei by using the A-core potentials described in Eqs. (pOj), (|3l|), ( |35|) and (36). 
In the case of the YNG and ORG potentials , we use the parameters given in Eqs. ([}(]) and (31). 
For the BOY potential we fit the parameter Vq (MeV) in Eq. (|36) whereas for the SW1 and SW2 
potentials, the best fit parameter is g aAA = 5o-aa/v / 47t in Eq. (3|). The errors (in brackets) in 
the fit parameters and in the mean squared radius are due to the experimental errors in the 
binding energies which are being fitted to. 
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Figure 1: A-core potentials (in MeV) for different hypernuclei. Solid lines correspond always to 
the BOY or SW1 potentials. For 9 A Be and ^ 09 Pb, the dotted lines stands for the SW2 potential, 
and for ^He we also give the YNG and ORG potentials. 
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SW1 


1 




2.4 


ORG 


SW1 


1 


3.25 ±0.03 


0.4 


YNG 


SW1 


1 


3.27 ±0.03 


2.7 


BOY 


BOY 


2 


2.37 ±0.02 


1.3 


ORG 


BOY 


2 


2.36 ±0.02 


0.03 


YNG 


BOY 


2 


2.37 ±0.02 


3.6 


SW2 


SW2 


2 


2.37 ±0.02 


2.1 


ORG 


SW2 


2 


2.36 ±0.02 


0.6 


YNG 


SW2 


2 


2.37 ±0.02 


4.1 



Table 6: HF results for g a \A obtained from a best fit to the -Baa data. The several rows 
correspond to different A-core potentials and different values of the cutoff AqAA- In all fits, the 
wAA coupling is kept fixed to g u AA = 29ujNn/3. The last column gives the x 2 by degree of 
freedom for each fit. Errors are statistical and have been obtained by increasing the value of % 2 
by one unit. 



determination of the A-core potential do not constitute an important obstacle to learn 
details about the AA interaction. 

Fits using short distance repulsive A- helium core potentials (YNG or SW1) provide 
significantly higher values of x 2 /dof than those obtained when the ORG or BOY poten- 
tials are used for helium. Besides, the only statistically significant contribution to the 
X 2 /dof comes from the helium datum. From this fact, one might conclude that short 
distance repulsive potentials for helium are not favored by the experimental value of the 
A \He binding energy. That contradicts the findings of Refs. |52L ^Pj where it was shown 
that short distance repulsive A-helium potentials were favored by the experimental value 
of the mesonic decay width of the A He hypernucleus. However it is difficult to draw any 



firm conclusion because, as we will see in Subsect. |4]|], YNG and SW1 potentials for 
helium are not so much statistically disfavored when r 12— correlations are included in the 
AA wave-function. The HE term in Eq. ([5]), which expected value vanishes for HF-type 
wave-functions but does not when r 12 — correlations are included, improves the simulta- 
neous description of an extremely light hypernucleus, as helium, and not as light ones, as 
beryllium or boron. 

For the sake of simplicity, when quoting our final HF results, we take for the central 
value and its statistical error the results obtained with a BOY A— core interaction for all 
hypernucleiQ. The results obtained with the rest of the potentials are used to estimate 
the size of the systematic error due to the A-core potential uncertainties. To be more 



3 The g<xAA values obtained when we use the ORG or BOY A-core potentials for helium are practically 
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specific, for a given value of the cutoff mass Ao-aa we take for this systematic error the 
statistical dispersion between the different values of g a AA obtained with each of the A— core 
potentials. Besides we also include another systematic error to account for the dynamical 
re-ordering effect in the nuclear core due to the presence of the second A, as we discussed 
in Sect. | after Eq. (g). 

Thus for the results of Table |6| corresponding to g^AA / QuiNN =2/3, we find 



/-_/ 3.23 ±0.03 ±0.02 ±0.03 A aA A = 1 GeV . , 

g*AA/ V47T - | 2 3? ± Q Q2 ± Q Q1 ± Q Q2 A ^ = 2 (54) 

where the first error is statistical and the second and third ones are the systematic errors 
due to our uncertainties in the A— core potential and in our treatment of the dynamics 
of the nuclear core, respectively. The latter one should account for an uncertainty of 
about 0.5 MeV, as discussed in Sect. in our theoretical determination of -Baa- For this 
systematic error we have taken an error of the same size as the statistical one. Thus, the 
inclusion of this error should account for an uncertainty of about 0.5 MeV (approximate 
size of the typical errors of the experimental data for ABaa) i n our determination of the 
total energies of the double— A hypernuclei. By adding all errors in quadratures, because 
all of them come from totally uncorrelated sources, we finally obtain 



, r— J 3.23 ±0.05 A (jA A = lGeV 

= 2.37 ± 0.03 AfjAA = 2 GeV (55) 



As mentioned above, in Ref. |^6[ the hyperon-nucleon scattering data were fitted with 
S^aa/s^aw =2/3, Ao-aa = 1 GeV and g cr AA/"\/47i : =2.138. This coupling constant is not 
directly comparable with that in Eq. (|53|), because the one determined in this work cor- 



responds to an effective interaction. We will give some more details in the Subsects. |475 
and 



The AA potential used up to here has some non-localities, namely the terms propor- 
tional to 

■•• "(VY(a,r)) L ■ V + Y(a,r)V 2 ] (56) 

in Eqs. (|24|) and (|25|) . Because of the small momentum of the A particles in the hypernuclei 
and because these terms are suppressed by powers of (m mes0 n/' r ''ibaryon) 2 , the contribution 
of these non-localities is small. In order to have a local potential to be plotted, we have 
dropped out these terms and redone the fits. We find fits of similar quality to those 
obtained with the non-local AA interaction and values of g a AA smaller than the previous 
ones by less than about 0.5%, indicating that the effect of the non-local terms is negligible. 

In Fig. H], we analyze the dependence of the AA interaction on the value of the ujAA 
coupling parameter, g^AA, an d the cutoff A ct aa- We show different local AA potentials, 



the same, and they only differ by less than one third of the statistical error of any of them. Thus, we 
have decided to quote those results obtained from a BOY A— helium core interaction in order to be able 
to use the same A— core potential for all hypernuclei studied in this paper. 
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obtained from best fits to the data. Despite the different shape and magnitude of the 
interactions shown in the figure, all of them give values for x 2 /dof of the order of 1, 
indicating the goodness of the fits and the impossibility of selecting any of them only by 
means of the binding energies, -Baa- The greater the parameter g w AA, the smaller the 
statistical and systematic errors in the fitted parameter g a AA are. Thus, the systematic 
(statistical) errors due to the uncertainty in the A-core potential vary from 2% (2%) for 
9uAa/9uiNN =1/3, to 0.2% (0.3%) for q u aa/ 9uNN =4/3. Also the x 2 /dof values decrease 
moderately for increasing values of the ratio guhkl 9uNN- 

As can be seen in Fig. ^] for a fixed value of the cutoff A ct aa, all potentials coincide at 
the same point. That implies^ that there is a linear correlation between the couplings 
^ctAA' &Iaa > defined in Eq. @, and (g£ AA - g£ AA ). Indeed we find excellent fits of the 
couplings to a dependence of the type 

9aAA = a + bx Caa + cx (^aa-Caa), (57) 

with approximately similar values for the parameters b and c, and certainly compatible 
within the statistical errors. Thus, there is an accidental cancelation of the dependence 
of Eq. (^) on g^ AA . To take advantage of this fact, we perform a two parameter fit and 
we find 

g 2 aAA _ i (2.97 ± 0.19) + (0.830 ± 0.011) x 
4vr \ (1.58 ±0.09) + (0.448 ±0.005) x 

for our local family of potentials. 

For the the full non-local AA potential, we 



o 2 

^AA 



X 2 /dof 
X 2 /dof 



o.oi, K 

0.02, K a 



AA 



AA 



lGeV 
2GeV 



(5* 



find 



^aa _ f (3.00 ±0.19) + (0.836 ±0.011) x x 2 /dof = 0.02, A ffA A = 1 GeV 

4tt [ (1.58 ±0.09) + (0.454 ±0.005) x ^aa, x 2 /dof = 0.02, A ctA a = 2 GeV 

(59) 

To perform the fits in both the local (Eq. (|58"D ) and non-local (Eq. fl59| )) cases, we add 
in quadrature, for each value of g^AA, the statistical and systematic errors of the fitted 
parameter g a AA- 

We would like also to mention that, here again, as we expected after our discussion in 
Eq. (P8"D, and in contrast to the findings of Ref. [4~5|, there is a strong correlation between 
the parameters g a AA and A ct aa and the quantity g<jAA x (1 — m2 / -^aa) remains constant 
within approximately 3%. This means that, within the HF frame of the analysis of AA 
hypernuclei, the use of different cutoffs A ct aa only amounts to redefine the g a AA coupling. 
Then the point at which the potentials coincide is approximately the same for both sets 
of AA potentials (A ct aa = 1 GeV and A ct aa = 2 GeV), as can be seen in Fig. |2|. 

10 Let us suppose that the point where all potentials coincide is ro and the common value of the potential 
is Vo- Then one finds 

rL = S=0' \ | T/i=S=0/ 



AA 



equation which automatically provides a linear relation between the couplings g^AAi S^aa an( l (SZ 
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Figure 2: Different 1 5o AA potentials, obtained from fits to -Baa- Left (right) figures correspond 
to results obtained within the HF (VAR) approach. In all figures the five curves correspond to 
the values 1/3, 1/2, 2/3, 1 and 4/3 for the ratio g^AA/duiNN- The cutoff mass A ct aa is 2 GeV 
(top) or 1 GeV (bottom). We have used a local AA potential, as described in the text after 
Eq. (|56|). BOY A-core potentials were used in all cases. 
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A ctA a = 1 GeV 

pert HF 
y erAA 


A ctA a = 2 GeV 

pert HF 

9<rAA-9 CT AA (%) r,f~f\ n KF 1 (07 r ) 

WctAA 


1/3 
1/2 
2/3 
1 

4/3 


5.2 1.8 
3.2 1.2 
1.9 0.9 
0.7 0.5 
0.3 0.3 


5.9 1.8 
3.7 1.2 
2.5 0.8 
1.0 0.5 
0.4 0.3 



Table 7: Comparison of perturbative (pert) and HF results for different types of AA interaction 
(the small non-local terms have been neglected). BOY A-core potentials have been used for all 
hypernuclei. The columns labelled by statfg^Aj stands for the statistical error (in percentage) 
of the fitted parameter g^ AA ■ 

4.3 A Perturbative Approach 

The fact that the quantity AB AA /2B AA is smaller than one for all hypernuclei considered 
in this work has made us to explore the possibility of using a perturbative scheme. Sur- 
prisingly, we find that in a good approximation the AA potential, Vaa, can be treated 
perturbatively, as the results of Table [7| indicate. In this table, we compare the HF-local 
results of Eq. (|58|) with those obtained from best fits to the data using the approximation 

AB AA = -(v AA (r 1 -r 2 )-^-^) , (60) 
\ M A / # , ) 



(0) 
AA 



where denotes the expected value of the operator O in the state $ and $ 

V 9 a( 1 "i) • ^a{^2)i where tp A is the A-wave function in the single-A hypernucleus A A 1 Z. To 
simplify we have also neglected the non-local terms in Vaa- Note also that the expected 
value of the HE term is zero in the state 3?a A . 

pert HF 

In the worst of the cases (g u AA/ ' 9uNN = 1/3 and Ao-aa = 2 GeV) the ratio < Lb \p aAA , 

9aAA 

which accounts for the difference between the HF and the perturbative approaches, is 
only 5.9% and its value decreases to only 0.3% for the case g^w/ QujNN = 4/3, being this 
latter value of the same order as the the statistical error (in percentage) of the fitted 
parameter g^ AA . Perturbative values of the fitted parameter g CT AA are always larger than 
the HF ones, though for most of the cases both sets of couplings are compatible within 
total (statistical and systematic) errors. The values of x 2 /d°f are about a factor of 1.5 
larger within the perturbative approximation than in the HF one. 

The fact that the perturbative and the HF schemes give the same results within 
approximately 5% can be used to understand the existence of strong linear correlations of 
the type of Eq. (|57|). For simplicity let us discuss only the case of the local AA potential. 
For each hypernucleus (helium, beryllium and boron), because Eq. ( |60D is satisfied in 



26 



a good approximation, one gets a linear relation between the couplings g^AA; qIjaa and 
S^AA' namely 



2 

9aAA 



-AB 



AA 



<K,V> 



(^T 5= 7^aa)*(o) <^T 5= 7^AaL(0) 



r 2 

3o;AA 



-*- A A 



where 



V 2 (r) 



(V l l = s =°/91aa)^ 



X (fi^AA — §Lya); 



V, 



L=S=0, 



r 



5wAA 



/(ft 



JojAA 



?uAa)' 



(61) 

(62) 
(63) 



and V^~ s 



a = a, cu) is the local part of V^ =5=0 . The fact that we find a reasonable 
simultaneous description of the three hypernuclei, with the same value of the parameter 
fi'o-AA for each value of ^AAj implies that the coefficients of the couplings in Eq. ( |6TD are 
rather independent of the nuclear core. 

Note that in the case of the non-local potential, linear relations between the couplings 
can be also obtained once the non-local terms of the AA potential are included in the 
right hand side of Eq. (^1|) . 



4.4 AA Interaction: Variational Results 

In Table ||, for given values of the ratio q^aa/ Q^nn and of the cutoff A ct aa, we show 
the values of the parameter g CT AA obtained from best fits to the Baa binding energies 
by using the variational method described in Subsect. 0. We use the same procedure to 
give the central values and the statistical and systematic errors as in the HF case. In 
all situations (i.e., different ratios g W AA/ g^NN, different values for the cut-off parameter 
Ao-aa and A-core potentials), we find statistically acceptable fits. We would like just to 
mention that the use of the ORG potential for helium provides smaller x 2 /dof than the 
BOY A-core potential used to quote the central value of the fitted parameter. Namely, 
values of x 2 /dof smaller than 0.1 for the ORG potential versus values around 0.5 for the 
BOY case. The HE term improves significantly the simultaneous description of all three 
hypernuclei when A-core potentials (BOY,YNG, SW1 or SW2) not as attractive at short 
distances as the ORG one are used for helium. To be more specific, for the case of the 
BOY (YNG) potential, if the HE term is not included the typical values for x 2 /d°f are 
of the order of 1.5 (4.0), between two or three times larger than those obtained when the 
HE piece is considered. In all cases the largest contribution to the x 2 /dof comes from 
the helium piece of data. 

Local and non-local potentials lead to different values for the fitted coupling ^aa- 
These changes vary from 1% to 7% when the ratio (Ju>aa / 9u>nn increases and they are 
significantly greater than the statistical errors (except for the lowest value of ^aa) and 
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S^aa/ 






A.aa (GeV) 


LOCAL 


NON-LOCAL 


1/3 


1 


1.89 ±0.03 ±0.02 


1.91 ±0.03 ±0.02 


1/2 


1 


2.33 ±0.03 ±0.02 


2.38 ±0.03 ±0.02 


2/3 


1 


2.76 ±0.03 ±0.02 


2.85 ±0.03 ±0.02 


1 


1 


3.55 ±0.03 ±0.02 


3.74 ±0.03 ±0.02 


4/3 


1 


4.26 ±0.03 ±0.02 


4.56 ±0.03 ±0.02 


1/3 


2 


1.397 ±0.018 ±0.008 


1.415 ±0.018 ±0.008 


1/2 


2 


1.764 ±0.017 ±0.007 


1.796 ±0.016 ±0.007 


2/3 


2 


2.145 ±0.016 ±0.007 


2.199 ±0.016 ±0.006 


1 


2 


2.876 ±0.017 ±0.006 


3.005 ±0.017 ±0.006 


4/3 


2 


3.551 ±0.019 ±0.006 


3.772 ±0.018 ±0.006 



Table 8: VAR results for g a \A obtained from best fits to the -Baa data. The several rows 
correspond to different values for the ratio guiAA/guNN and different values of the cutoff A ct aa- 
Central values and their statistical uncertainties (first set of errors) are obtained using BOY A- 
core potentials for all hypernuclei. The x 2 /dof values for all cases are about 0.4-0.5. Statistical 
errors are obtained by increasing the value of x 2 by one unit. The second set of errors account for 
the dispersion in the results due to the use of different A-core potentials, as it was discussed in 
Subsect. W3. To account for dynamical re-ordering effects in the nuclear core due to the presence 



of the second A, further systematic errors should be added. As it was discussed in Subsect. 4.2 , 
the statistical errors provide a reasonable estimate of these additional systematic errors. We 
give results for the full AA interaction (non-local) and the local reduction of it discussed in the 
text after Eq. (56). 
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much more appreciable than those found within the HF scheme. Variations of the cutoff 
mass A ct aa lead to changes in the quantity g a AA x (1 — ^/A^aa) fr° m 3% to 10% when 
the ratio g^AA/ g^NN increases from 1/3 to 4/3. These changes are again much more 
pronounced here than in the HF case, indicating that the variational wave functions 
for the two A's system lead to bigger values of the transfered momentum between both 
particles. This is also consistent with the fact, mentioned above, that within the VAR 
scheme non-local and local potentials lead to differences much more significant than in 
the HF one. 

In Fig. |] (right) we show the local AA potentials corresponding to the different com- 
binations of the ratio g^AA/g^NN and the cutoff mass parameter A ct aa given in Table ||. 
In this figure, for comparison, and in the same scale, the HF local potentials discussed in 
the previous section are also shown. VAR potentials are significantly less attractive than 
HF ones, as the variational principle ensures. The larger the ratio #u>aa/<?u>aWj the bigger 
this effect is. 

For both local and non-local VAR potentials, we find that the coupling constants are 
correlated in the form given in Eq. Q57|), but with poorer x 2 Jdof values (with x 2 /dof 
ranging from 0.12, for a local AA interaction and A ct aa — 1 GeV, to 0.8, for a local AA 
interaction and A ct aa = 2 GeV) than in the HF case. But, for this variational scenario, the 
values of parameters b and c of Eq. (|57| ) are not similar, and a linear relationship between 
g^AA and 5^aa can n °t be casted. Therefore, unlike the HF families, the obtained families 
of variational potentials do not coincide at a single pointQ. 

The existence of a relation of the type of Eq. ( p7|) also in the VAR scheme can be un- 
derstood if one applies the same type of considerations used in the discussion of Eqs. ([|]) - 
TTf) to the relation 

Baa = - ( E f + VaaCtS) ] + V^fi - 3) - . (64) 

\i=l,2 V 2 t*A J M A I R 

AA 



In Fig. [3] we show graphically the differences between the fitted parameter g a AA ob- 
tained in the HF and VAR schemes for different values of g^AAj different cutoffs A ct aa 
and using local (dotted lines) and non-local (solid lines) AA interactions. The non-local 
part of the AA potential is repulsive, as can be deduced from the figure. However, the 
differences between VAR local and non-local parameters are always much smaller than 
the differences between HF and VAR parameters. 

Once we have seen that within the HF approach (or even in the perturbative approx- 
imation) we could described reasonably well the experimental data, one might question 
about the need of performing a variational study of the double— A hypernuclei. To check 
this in Table we compare, for a fixed AA interaction, the results of the perturbative, HF 

11 The equation dV\A(r) / dg^AA = has as solution a value of r independent of g^AA (what implies 
that all potentials coincide in a single point) only if there exist a linear relation between the couplings 
9aAA an d 9wAA an d a small term proportional to {m^ /ttia) 2 g^AAfujAA, in the to— exchange part of the 
potential is neglected. 
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Figure 3: Different couplings g a kk/v^ versus g^w / g^ nn found for the Sq AA potentials 
obtained from best fits to .Baa data within the HF and VAR approaches. The cutoff mass A ct aa 
is 1 GeV (left) or 2 GeV (right). In both plots the solid (dotted) lines correspond to a non-local 
(local, as described in the text after Eq. (|56|)) AA interaction and the five points with errors 
correspond to the values 1/3, 1/2, 2/3, 1 and 4/3 for the ratio guiA\/guNN- The lines have been 
drawn only for guiding the eye. BOY A-core potentials were used in all cases. For comparison, 
the square in the left figure stands for the potential of model A of Ref . [ 46 ] . 



and VAR approaches for the binding energy of We also show the result obtained 

if no correlations (i.e., dependence on ri 2 ) are included in the variational wave function 
(that corresponds to fix c = in the variational basis of Eq.(|10|)). Even in this latter 
case the VAR scheme does not reduce to the HF one because the wave function can not 
be gathered as a product of a function of r\ times the same function now of r 2 . The 
HF approach can be considered as a variational one, where the space of functions are 
limited to those which can be written in a factorizable way. The numbers of the ta- 
ble show that enlarging this space to include some non-factorizable wave-functions (but 
without dependence on r 12 yet) does not modify drastically the HF results. However, the 
inclusion of r 12 — correlations has a drastic effect in the binding energy. This can be also 
appreciated in Fig. [|, where the dependence on r 12 of the AA potential (used to obtain 
the numbers of Table |9|) and the HF and VAR densities of probability of finding the two 
A's at a relative distance r i2 (P(r i2 ), denned in Eq. (0)) are shown. The effect of the 
inclusion of ri 2 — correlations is twofold. First, it reduces the probability of finding the 
two particles at very short distances, where the potential is highly repulsive, with values 
at rq 2 = as large as tens of thousands of MeV. Second, it increases the probability of 
finding the particles at relative distances of the order of 1 fm, where the potentials have 
their maximum attraction. Thus, as a net effect an important reduction of the expected 
value of the energy of the system is provided. 

The AA interactions obtained by a — uj exchange behave almost like a hard-core at 
short distances, then in a very small region and close to the origin (ri 2 « 0.5 fm), they 
pass from being extremely repulsive to reach their maximum attraction. On the other 
hand, note that due to the phase space P(r i2 ) is suppressed at small values of 7q 2 by 
at least a factor r\ 2 . Thus the HF two body probability P(ri 2 ) takes very small values 
in the region where the potential is repulsive but it also misses the maximum attraction 
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B AA [f A B] [MeV] 



PERT. 27.3 

HF 27.8 

VAR NO CORR. 28.4 

VAR 36.9 



Table 9: Comparison of the results obtained for -Baa i n boron double— A hypernuclei, using 
different approximations. In all cases the same local AA interaction (with parameters g a AA/V^, 
9luAa/9ujNN and A^aa fixed to 3.22, 2/3 and 1 GeV respectively) and BOY A-core potential were 
used. The coupling g<jAA wa s obtained from a best fit to the data within the HF approach (see 
results of Eq. (58)). The first, second and fourth rows correspond to the standard perturbative, 
HF and VAR calculations. The third (labeled VAR NO CORR.) row corresponds to a variational 
calculation where the ansatz wave function does not depend on r\2 (that corresponds to fix 
c = in the variational basis of Eq.(pX|)). The experimental binding energy for boron double— A 
hypernuclei is 27.5 ± 0.7 MeV. 



/ 1 2 3 4 5 6__ 

Vi 0.9727 0.02280 0.00200 0.00075 0.00054 0.00040 0.00028 
E l k = Vk 0.9727 0.9955 0.9975 0.9982 0.9988 0.9991 0.9994 



Table 10: Probabilities Vi (defined in Eq. (|i~2|)) for several waves. Results were obtained for 
boron double— A hypernucleus. A local AA interaction (with parameters go-AA/V^r = 2.76, 
9u>aa/ QujNN = 2/3 and AoAA = 1 GeV, see Table |8|) and BOY A-core potentials were used. The 
errors are always less than one unit in the last digit. 

region of the AA interaction. The inclusion of a dependence on r 12 in the variational 
wave function allows both to keep small (and even smaller than in the HF case) the two 
body probability in the hard-core region and to increase this two body probability in 
the adjacent region where the potential reaches its maximum attraction. The steeper 
AA potential, the bigger the difference between the VAR and HF predictions becomes. 
Thus, the effects of including r\i— correlations are much more pronounced for values of 
the ratio q^aaI 9^nn bigger than the standard 2/3 one used in this discussion. Therefore, 
the VAR approach, including correlations, provides the same binding energies than the 
HF approach with significantly less attractive AA interactions, as we showed in Figs. |2] 
and |3|. 

To finish this section we would like to devote a few words to the multipolar content 
of the VAR wave functions. In Table |1(] we show the probability Vi, defined in Eq. (|12]) , 
of finding each of the two A's with angular momentum I and coupled to L = for 
boron. Although a particular AA interaction is used, the gross features of the multipolar 
decomposition do not depend significantly on the selected interaction. Only few multipoles 
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> 




Figure 4: Top: HF (solid) and VAR (dotted) probabilities of finding the two A's at a relative 
distance r\z as a function of this distance. Computations were done in ^a.B and using the AA 
interaction (bottom figure) specified in Table |9[ 
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n 1 

9aAA/ 


/ A ~~ 
V 47T 


^o-AA 


/ 1 A—. 

j V 47T 






! A 

\P~trhK — 


1 CjeV) 


/ A 

l^crAA " 


- 2 GeV) 




A MA [GeV] 


HF 


VAR 


HF 


VAR 


without 




3.23 ±0.03 


2.85 ±0.03 


2.37 ±0.02 


2.199 ±0.016 




1.5 


3.38 ±0.03 


2.85 ±0.03 


2.48 ±0.02 


2.232 ±0.018 


with 


2.0 


3.52 ±0.03 


2.74 ±0.03 


2.58 ±0.02 


2.16 ±0.02 




2.5 


3.60 ±0.03 


2.62 ±0.03 


2.64 ±0.02 


2.06 ±0.02 



Table 11: Best fit results for the uAA-coupling obtained from different non-local 
AA— interactions for both HF and VAR approximations. We compare results for both a + lo 
(already presented in Tabs. |6| and |8| ) and a + lo + 0— interactions. BOY A-core potentials were 
used for the three double— A hypernuclei. The parameter g u AA/9u>NN has been fixed to 2/3. 
Furthermore, the 0— couplings have been fixed, by means of S'C/(6)-symmetry, (see Eq. (|65|)). 
Errors are only statistical. 



contribute^ being / = the dominant one and contributing appreciably only up to the 
waves I = 2 or I = 3. The higher the orbital angular momentum of each A, the farer from 
the origin the corresponding wave function is and therefore the smaller the overlap of the 
two-body wave function with the attractive part of the potential is. 



Results of Table ITUl which show predominance of the I = component in the VAR 



wave-function, do not contradict the fact that correlations are important to lower the 
energy of the double— A hypernuclei. Indeed, the hamiltonian is not diagonal in the basis 
with well defined single-particle orbital angular momentum. Hence, non-diagonal matrix 
elements contribute significantly to the hamiltonian expectation value. 



4.5 Contribution of the 0— Exchange. 

In this Subsection we discuss qualitatively the effect of the 0— exchange in the AA potential 
in the medium. The 0— potential can be obtained from that of the lo— exchange given in 
Eqs ( P5D - fl2"?D by the obvious substitutions: m u m^, and [<?,/, A]^aa — > [g, f, A]^aa- 
In the spirit of the Bonn-Julich models we use SU(6) symmetry to fix the 0— couplings, 

fi^AA = —9cjAa/ V2 

UaA = V^fwAA- (65) 



As it was discussed in Subsect. |2.1| because within this model the meson does not couple 
to nucleons, we do not have much information about the cutoff-mass A^aa. Assuming 
that this cutoff should be similar to that entering in the a;AA— coupling and certainly 



12 Note that we show probabilities: the components of the wave function for each value of I are given, 
up to a sign, by the squared root of the numbers presented in the table. 
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bigger than the meson mass, we have studied three different values, 1.5, 2 and 2.5 GeV. 
Results are presented in Tables |ll] and [12|. We would like to make a few remarks, 



The effect of the 0— exchange depends strongly on the value of the cutoff mass A^aa, 
such that this effect passes from being small (A^aa = 1.5 GeV) to be significant 
(A^aa = 2.5 GeV). Indeed, this is easily understood by looking at the top graph of 
Fig. |5] where the uj and the three uj + 4> potentials considered in Table [II] are plotted 
(neglecting the small non-local pieces). The several uj + potentials significantly 
differ in the region of relative AA distances ranging from 0.1 to 1.3 fm, where the 
product of the AA-interaction times the probability (P(ri2)) of finding the two A 
particles at relative distance Tyi takes non negligible values. 



As the numbers of Table [12] indicate, the 0— contribution is smaller than the uj me- 
son one and at maximum it might be of the order of 20% of the latter one. This 
is not surprising, because as we mentioned earlier, the exchange of heavier mesons 
than the a and u ones, essentially modify the short range behavior of the inter- 
action, and hence becomes more significant for processes where larger momentum 
transfers than those involved here are relevant. On the other hand, by the use of 
different values of the a— cutoff and also different values of g^AA we have already 
explored a wide range of short distance behaviors of the effective AA— interaction 
in the medium. The above discussion, together with the very scarce set of data 
available, have prevented us to include the 0— meson contribution in the present 
determination of the effective AA interaction. Its inclusion would require, even as- 
suming SU(6) — symmetry, to deal with at least one more free parameter, and it 
might obscure the analysis presented in this work (differences between HF, pertur- 
bative and VAR approaches, size of systematic errors ■ • ■). On the other hand, the 
possible effect of the 0— exchange potential in the tail of AA wave-function, though 
small, is mostly reabsorbed by the use of effective a + uj— couplings and thus, the 
inclusion of the contribution in the potential would not modify significantly the 
mesonic decay of the double— A hypernuclei studied in the next Subsection. How- 
ever, when trying to simultaneously analyze the double— A hypernuclei data and the 
S = —2 baryon dynamics in the free space, one should not neglect this contribution 
because it might lead to sizeable changes in the values of the o— coupling, as the 
numbers presented in Table [TT] clearly show. 

A distinctive feature of the results presented in Tables [TT] and [T2|, is that the char- 
acter of the 0— contribution though attractive within the VAR scheme, is repulsive 
within the HF one. Hence, the differences between the HF and VAR approaches get 
amplified. 0— exchange does not lead to a purely repulsive potential, as it is the case 
for the uj— meson, because the ratio f / g is positive, instead of negative, and bigger 
(in absolute value) for the 0— potential than for the uj— one. Thus, the 0— potential 
though much more repulsive than the uj— one at short distances, becomes however 
attractive at distances bigger than let us say 0.4 fm (see Fig. |^ Q). In the case of 
the HF approach, though the probability P(ri 2 ) is quite small at short distances, 



13 Note that in the limit A^aa ~* oo, the repulsive core becomes a S peak at the origin. 
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A MA [GeV] 


HF 


VAR 


1.5 


1.95 


-0.002 


2.0 


3.68 


-1.14 


2.5 


4.53 


-2.33 



Table 12: Expected values of the local part of the 0— potential, for several (f>AA cutoffs, in the 
AA state corresponding to the HF and VAR solutions (with no 4> contribution and A ct aa = 1 
GeV) given in the first row of Table |ll]. 4>— couplings have been fixed by means of SU(6)- 
symmetry. For comparison, the expected values of \ V^\ u+w) are 14.68 and 10.02 MeV for the 
HF and VAR wave-functions respectively. 

because the huge and repulsive values taken by the 0— potential close to zero, the 
balance between the repulsive and attractive contributions favors the former ones. 
The inclusion of r 12 — correlations in the VAR scheme allows to have simultaneously 
smaller and bigger values of P(n.2) than in the HF case in the regions where the 
<p— potential is repulsive and attractive respectively. As a net effect, the balance 
between the repulsive and attractive contributions favors now the latter ones. This 
is easy to appreciate in Fig. ^[ 



4.6 Mesonic Decay and Binding Energies of Double— A Hyper- 
nuclei 

One might think that the mesonic decay of these double— A hypernuclei could depend 
significantly on the details of the effective AA interaction and thus it could be used to 
differentiate between the different potentials shown in Fig. |2|. There are some theoretical 
uncertainties^ in the calculation of the mesonic decay of A He ^ |29|, gl], related to the 
nuclear core-A interaction, and this has prevented us from looking at \\B.e to check the 
dependence of the mesonic decay on the details of the AA interaction^. Thus, we have 
looked at the case of aaB- F° r this hypernucleus, we find that the pionic decay width of 
12 A B changes only around a 5% when different A-core potentials (BOY, SW1, SW2) are 
considered. 



The mesonic decay has been computed following the method exposed in Sect [| We 
use the ten VAR non-local potentials presented in Table BL and we find that the mesonic 



14 As can be seen in Fig. [T] and in Table [|, different A-core potentials for helium lead to significantly 
different mean squared radius of the A orbiting around the nuclear core. The mesonic decay process is 
quite sensitive to the tail of the A wave-function and then it depends strongly on the details of the A-core 
potential used. 

15 In Refs. |2^] and |39| the mesonic decay of the double— A hypernucleus A A^- e nas b een calculated. 
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Figure 5: Top: uj— (solid line) and u + 4>— local pieces of the AA potentials predicted by 
S'L r (6)-symmetry. Three different values (1.5, 2.0 and 2.5 GeV) for the cutoff A^aa have been 
used. The larger the cfiAA cutoff, the bigger is the resulting repulsion at the origin. Bottom: 
Two different set of curves are being plotted as a function of the AA relative distance, r\i : a) 
HF and VAR probabilities, obtained from the potentials specified in the first row ("without </>" 
and A CT AA = 1 GeV) of Table 11, of finding the two A's at a relative distance r\i , b) the product 
of these probabilities times the A^aa = 2.0 GeV local (^—potential used in Table This latter 
set of curves have been divided by a factor 50 and the area below these curves gives, up to a 



factor 50, the expected values given also in Table 12 



36 



decay width varies by 5% at most, making then this quantity unappropriated to choose 
between the different potentials discussed above. However, this fact allows us to predict 
the mesonic decay and, by using the potential of Eq. ( fT/) ) in addition to Cohen-Kurath 
spectroscopic factors for the lp— shell to describe the n B, 12 C* and 12 B* nuclear cores, we 
find 



m 3 A B->X + A + 7r°) 



= 0.062 ±0.002 ±0.002, (66) 

J- A 

r(l 3 A B -> X + A + 7T-) 

- = 0.270 ±0.008 ±0.008, (67) 

where Ta = ± T^l e is the total decay width of the A in the vacuum, with T^} ee 
defined in Eq. (pffi). The central values are obtained with ^aa/ QluNN = 2/3 and Ao-aa = 1 
GeV. The first set of errors accounts for the statistical error of g a AA whereas the second 
one accounts for the spreading of the results obtained with the different combinations of 
ratios c/uiaa/ QloNN (1/3, 1/2, 2/3, 1, 4/3) and values (1 and 2 GeV) for the cutoff A ct aa- 



Further systematic errors should be added to those quoted in Eqs. (pq ) - (|67| ) due to 
both the uncertainties in the A-core potential (error of about 5%, as mentioned above) 
and the no inclusion of the nuclear core distortion effects (errors of about of the same size 



than the statistical ones, as discussed in Subsects. [O] and |4.4j ). Adding all the systematic 
and statistical errors in quadratures we end up with total errors of about 7% for both the 
7T° and 7r~ decay widths. 

The main (> 95%) contribution to these decay widths comes from processes where the 
outgoing nucleon is in a bound state of the daughter hypernucleus, processes described in 
Eqs. © - ©. 

Note that the final nuclear state is described by means of a simple shell-model sup- 
plemented by an effective interaction (Cohen-Kurath) for the lp— shell nucleons. Results 
presented above might depend on the specific details of the used central and residual inter- 
actions. This is a problem, which is not specific of the mesonic decay of double— A hyper- 
nuclei and it is already present in the studies of the mesonic decay of single A— hyper nuclei. 
Indeed, if the residual interaction is neglected, we obtain with the potential of Eq. (f4?j) 
values for decay widths (in units of T A ) of 0.074 ±0.004 ± 0.004 and 0.370 ± 0.008 ±0.008 
for 7T° and tt~ decays. The difference between these values and those quoted in Eqs. (^) - 
( |6TD might give us an estimate of the size of our uncertainties. A recent calculation of the 
7r~ decay [fiO" , where the nuclear states are been constructed in the shell model within 



p— shell configurations with Cohen-Kurath interaction [p7| , gives 0.325. This value fits 
within our range of results. 

To minimize the effects of this new source of systematic errors, it is interesting to 
define the ratios 
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ra 3 A B^x + A + 7r ) 



2r( 1 2B -> X + 7T°) 

r( A 3 A B -> X + A + 7T- 
2r( 1 iB^X + 7r-) 



0.63 ±0.02 ±0.02 ±0.02, (68) 
0.67 ±0.02 ±0.02 ±0.02, (69) 



which differ from the naively expected value 1. These values have been obtained by using 
the potential of Eq. fl47f) and Cohen-Kurath Ip— shell spectroscopic factors to describe 
the nuclear cores involved. Neglecting the residual interaction the potential of Eq. (f47j) 
leads to similar values (0.65 and 0.72 for tt° and n~ decays respectively) for these ratios, 
as we expected. The meaning of the first two set of errors in the equations above is the 
same as in Eqs. (|66"D - (|67|) . The last set of error accounts for the systematics due to the 
no inclusion of the nuclear core distortion effects. On the other hand, in these ratios, the 
uncertainties due to the details of the A-core potential mostly cancel out. 

The mesonic decay depends both on the momentum carried by the outgoing nucleon 
(the greater the nucleon momentum, the less effective the Pauli suppression is) and on 
the tail of the overlap function defined in Eq. (|5T|). The first factor is totally determined 
by the energy balance in the reaction, which is the same for all potentials used. Thus, the 
fact that the mesonic widths calculated in the VAR scheme do not depend appreciably 
on the specific potential means that all overlap functions have a similar behavior at large 
distances. This is shown for two potentials in the bottom plot of Fig. ^. 

On the other hand, HF wave functions lead to results similar, within a 5%, to those 
presented in Eqs. ( |66D - (^) and one would again expect that both HF and VAR overlap 
functions have similar tails. This can be seen in the top plot of Fig. [| Finally, in this 
figure we also compare the HF and VAR overlap functions with the A wave function 
in the daughter hypernucleus (v?a)- Because the ratio -Baa/2-Ba is greater than one, 
the exponential decay at large distances of <£a is less pronounced. This, together with 
the change in the energy balance of the reaction, which makes more effective the Pauli 
blocking for double- than for single-A hypernuclei, explains why the ratios defined in 
Eqs (|j8|) - (|69"D are smaller than one. 

To finish this section, in Fig. |7] we show the A dependence of the ratios A.Baa/2-Ba and 
T(fX 2 Z ->• X + A + 7T '-) /2r( A A f X Z -»• X + 7T '-). As can be seen in the figure, we are able 
to give accurate predictions for both ratios which also turn out to be rather independent 
of the details, g u AA and A^aa, of the effective AA interaction. 

In Ref. [38| the A-dependece of ABaa has been also studied within a Skyrme-Hartree- 



Fock approach. There, it is also found that ABaa is a decreasing function of A. We find 
values for A£?aa within the broad range of possible results quoted in that reference. 



4.7 Nuclear Medium and Free Space AA Interactions. 

It is essential to state that the AA potentials determined in this work, effectively de- 
scribe the dynamics of the AA pair in the nuclear medium, but they do not describe 



38 





Figure 6: The projection of the two-particle AA state in a double— A hypernucleus, 3>aa> over 
the one-particle A state, ip\, in the corresponding single-A hypernucleus is shown versus r for 
Boron. This projection is the overlap defined in Eq. ([H]). In all cases, the A-core interaction 
is of the BOY type and the cutoff parameter A ct aa is 1 GeV. The different lines correspond to 
different calculations of <1>aa as follows. Top: Perturbative (solid), non-local HF (dotted) and 
VAR (dot-dashed) results, with g w AA./ 9uNN = 4/3 for the two last ones. Bottom: Perturbative 
(solid) and non-local VAR with two values for the ratio guAk/gwNN-, 1/3 (dotted) and 4/3 (dot- 
dashed). In the VAR and HF approaches g a AA is obtained from the overall fit to the -Baa data. 
Within the perturbative scheme, the overlap (^(^I^AA^ ^2)) coincides with <PA(r) (A wave 
function in the { A+ \Z) hypernucleus). 
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computed for different double— A hypernuclei. In the case of the ratio of mesonic widths with a 
7T° in the final state, the results have been multiplied by a factor 1 /2. BOY A-core potentials have 
been used for all hypernuclei and the non-local AA interaction associated to g^AA = ^9wNn/^ 
and A ff AA = 1 GeV has been used to determine the central value (see Table ||). Statistical and 
systematic errors in the fitted parameter g a AA lead to the error bars in the shown ratios, as it is 
explained in the text after Eqs. (^) - (^). Experimental errors in Ba have not been considered. 
In some cases the size of the errors (systematic and statistical errors are added in quadrature) 
is smaller than the symbols. 
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their dynamics in the free space. Indeed, what has been determined in this work is an 
effective interaction in the medium. This effective interaction is usually approximated by 
an induced interaction |5^, [75], |76[ (V^aO which is built up in terms of the A A — > A A 
(G A a), AN -> AN (G AN ) and AfJV -> NN (G NN ) G-matrices, as depicted in Fig. |. 
The induced interaction, Vfff, has the virtue of combining the dynamics at short dis- 
tances (accounted by the effective interaction Gaa) an d at long distances (polarization 
phenomenon) which is taken care of by means of the iteration of the particle-hole (ph) 
excitations (RPA series) through the effective interactions Gaa and GViv- The Gaa— , 
G\n— and the Gat at— matrices can be obtained from the S = —2, S — — 1 and S = 



baryon-baryon interactions in the free space |17], [24], [75]. In the 5* = —2 channel and 
near threshold, as discussed in the introduction, one needs to solve the AA — EN coupled 
channel G-matrix equations^} In a nuclear medium and with total energies of the order 
of 2mA, the AA — EN coupling might be suppresed due to Pauli blocking and thus the 
data of double— A hypernuclei would probe primarily the AA diagonal element^, V[ r K ee , 
of the A A — EN potential. Thus, Gaa might be roughly approximated by V[™ e . 



To establish the connection between V£jf and V[™ e is the aim of a future work p5 



and it is out of the scope of this work. As discussed above, it requires to understand the 
renormalization of the a, u propagation (it will be necessary to take into account that 



both carriers can excite ph components on their propagation through the nucleus |57| , |68fl ) 
and the role played by the EN intermediate states in the nuclear medium. Note that, 
within the Bonn-Julich model, the 0— meson does not couple to nucleons, and hence its 
propagation in the nuclear medium is not renormalized. 

In Table [T3| we present the smallest values of the g CT AA parameter which lead to 1 S Q A A 
bound states in the free space, when a + u + <p— AA interactions, with different cutoffs for 
the (f> potential, are considered. For the sake of simplicity, we have fixed q^kk/ g^NN = § 
and A ct aa = 1 GeV. (ft— couplings are given in Eq. (|65|). As can be deduced from the 



numbers quoted in this table and those already presented in Table |TT] the derived HF AA 
potentials would lead to AA bound states. However within the variational scheme in all 
cases the potentials are not attractive enough to bind the two hyperons0. 



Preliminary results of Ref. J55 show that the AA interaction in the medium, V^jf, is 
more attractive than that in the free space, V[ r ^ ee . Variational results (our best results) 
provide effective interactions unable to bind the AA pair in the 1 So channel. Therefore we 
can discard the existence of AA bound states in the free space when only the AA diagonal 
element, V/^ ee , of the AA — EN potential is considered. An independent confirmation 
of the no existence of AA bound states is given by the a— coupling (^oAa/V^t = 2.138) 
extracted from the study of hyperon-nucleon scattering processes which is well below of 



any of the a— couplings quoted in Table [13]. However, the AA — EN coupling, though likely 

16 The AA — EN mass difference is about 25 MeV. The mass difference between the AA and EE pairs 
is about 150 MeV. Thus one can safely neglect the influence of the latter channel in the study of the AA 
dynamics at threshold. 

17 It is to say that the contribution of the AA — > EN — > AA transition would be much smaller than 
the direct AA — > AA one (with no EN intermediate states and accounted for by the V/^ 66 piece of the 
S = — 2 potential), in this context. 

18 The conclusion is the same for the case A ct aa = 2 GeV. 
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without <p 


with (f) 


A MA [GeV] 




1.5 


2.0 2.5 


S^aa/ V^tt 


3.098 


3.105 


3.003 2.889 



Table 13: Smallest values of the g a AA— coupling leading to 1 So AA bound states in the vacuum. 
Different non-local interactions, all of them with g^Kh/ 9ujNN = § and A ct aa = 1 GeV, without 
and with the inclusion of the 4>— meson exchange potential, which coupling to the hyperons is 
fixed by means of the SU(6)— symmetry, have been considered. In the latter case three different 
cutoffs have been studied. 

suppresed in the nuclear medium, might contribute significantly in the free space. In the 



model proposed in Ref . |37] , the inclusion of the above coupling leads to more attractive 
interactions in the S = —2 sector, allowing for values of the AA 1 Sq free space scattering 
length as large as that of the nn system. Thus, to draw any firm conclusion about the 
existence of bound states in the S = — 2 baryon-baryon sector, a combined studied of the 
AA and EN systems is needed. 

We would like to stress that the important issue now is to clarify whether or not 
the difference between the g a AA/V^n values of 2.138 (obtained from hyperon-nucleon 
scattering data in Ref. P6[ ) and of 2.85 (obtained in this work for <7 w aa/ ' 9uNN — 2/3 
and Ao-aa — 1 GeV within the VAR scheme^ can be explained in terms of medium 
renormalization effects and/or the contribution of EN intermediate states to the induced 
interaction V^ A d and/or the contribution of heavier mesons not included here, such as 



the 4> meson. As discussed in Subsect. [4~5| , the inclusion of the 0— exchange piece of the 
AA potential reduces in some extent the difference between the g a AA couplings mentioned 
above. For instance, if a A^aa cutoff of 2.5 GeV is used, one finds values for ^oaa/v^tt of 
the order of 2.6 (Table 0). 



Finally a word of caution must be said here. The Jiilich group has also shown |77|, |7H 
that the correlated 2n and KK exchanges lead to a value for the ratio g a AA/gaNN of 
0.49. This value is much smaller than the results obtained in any of the present YN 
models E3, E3, |45|, EBl WA, ESI, for which this ratio varies from 0.58 to 1. Therefore, the 



intermediate range attraction in the AA?" and AA channels will be significantly reduced if 
one assumes the results of the microscopical calculations of Refs. [|77|, [78| . Then, a value 
for the ratio g u AA/gujNN smaller than 2/3, predicted by SU(3) and commonly accepted, 
would be needed to reproduce the AN scattering data. Thus, once all nuclear medium 
effects are understood, it will also be worth studying if a new set of parameters with 
guAA/gaNN = 0.49 and g w AA/ g^NN < 2/3 provides a simultaneous acceptable description 
of both double— A hypernuclei and Ap scattering data. 



1!) 



These values for the ratio g^AA/ 1 gu>NN and cutoff A ff AA correspond to those used in Ref. 
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Figure 8: Feynman diagrams relating the AA — > AA,AN — > AN and NN — » NN G-matrices 
to the induced interaction (VJ^jf) in the nuclear medium, which approximately accounts for the 
effective interaction determined in this work. 

5 Conclusions 

The dynamics of the AA system can not be studied directly in the free space. A p scattering 
provides an indirect source of information about it. Thus, the physics of double— A 
hypernuclei provides an independent and valuable method to confirm the findings obtained 
in the studies of hyperon-nucleon scattering. From these studies, it is commonly accepted 
that there exist three kinds of realistic models for the AA diagonal element of the AA — HiV 
potential, those developped by the Nijmegen, Tubingen and Julich groups. The Jiilich 
OBE potential for the hyperon-hyperon interaction is a natural extension of the Bonn 
model for the NN- interaction In this paper we have used for the first time Julich 
type AA OBE potentials to study the dynamics of double— A hypernuclei. We have found: 



Effective AA interactions in the medium extracted from Julich type models provide 
a simultaneous fairly good description of the binding energies of the three known 
double— A hypernuclei. 

Uncertainties in the exact nature of the A-nuclear core dynamics lead to uncertain- 
ties in the determination of the effective AA potential in the medium smaller than 
those due to the experimental errors in the fitted quantities ASaa- Thus, the ambi- 
guities in the determination of the A-core potential do not constitute an important 
obstacle to learn details about the AA interaction. 



The HF and the perturbative approaches, discussed in Subsects. fl^ and fO| respec- 
tively, lead to quite similar effective AA potentials. 

The inclusion of AA correlations in the variational scheme provides a better un- 
derstanding of the dynamical features of the system. It also enables less attractive 
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effective AA potentials than in the HF or perturbative approaches and more similar 
to those extracted from hyperon-nucleon scattering to describe the binding energies 
of the double— A hypernuclei. 

• Only few waves (up to I = 3) contribute appreciably to the multipolar expansion of 
the VAR wave functions. 

• The inclusion of the 0— exchange might be relevant to understand the effective AA 
interaction derived in this work in terms of the free space one. 

• The existing double— A hypernuclei data can not conclusively favour any particular 
choice neither of the ratio g^kh/g^NN in the interval [1/3,4/3] around the SU(3) 
prediction 2/3, nor of the cuttoff A ct aa in the range = 1-2 GeV. Thus, we end up 
with a whole family of effective AA potentials describing the ground state binding 
energy of the three known double— A hypernuclei. 

• The mesonic decay of \\B and the binding energies and pionic decay widths of heav- 
ier double— A hypernuclei (not discovered yet) turned out to be rather independent 
of the details of the effective AA interaction within the family of potentials described 
in the previous point. This fact has allowed us to predict accurately, for the very 
first time, the mesonic decay widths of medium and heavy double— A hypernuclei. 

• We have discarded the existence of 1 S'o A A bound states in the free space if the 
AA — EN coupling is negligible. 

The natural continuation of this work [|5lJ is the study of the nuclear medium mod- 
ifications of the AA interaction. Such study, together with the effective AA interactions 
derived here, will allow us to better understand the dynamical features of the AA interac- 
tion, at low energies, in the vacuum. Thus, it will be possible to contrast features of such 
interaction from two different and independent sources of data: hyperon-nucleon scatter- 
ing and AA hypernuclei. Statisticaly improved and new experimental data on double— A 
hypernuclei, exploring not only light but also medium and heavy nuclei, will be extremely 
valuable to achive such an objective. 
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Appendix 

Matrix Elements in the VAR Scheme 



The basic integrals to compute the matrix elements of the hamiltonian of Eq. (]5|) in the 
basis defined in Eq. ( |i~0|) are: 



Ti mn (a,c) 



■57T 



+00 



+00 



dr 2 r 2 m 



ri+r 2 



ri— n\ 



(70) 



with a, c and /, m, n real positive and integer numbers respectively. To compute the matrix 
elements of the A kinetic terms, of the HE term, of the A-core potentials and of the local 
part of the AA interaction, only non-negative values of l,m and n are needed. In this 
case, all the integrals can be obtained from Tooo by means of a recursion formula proposed 



in Ref. 55 



Blmn 

with r 000 given by: 



2^ r ^ l—l,m,n ^ r^m— l,n -^imn] 



a + c 



(47r) 2 (m + n)\ 
(a + c ) m+n+1 



ooo 



(4tt) S 



2a(a + c) 2 ' 



(71) 
(72) 

(73) 



To compute the matrix elements of the non-local part of the AA interaction, integrals of 
the type r^ m) _i, with l,m = 0,1,2 - •• are also needed. These can be obtained from a 
recursion as well fl66| , 



^l,m,-l — 

with T o -i given by 



1 

2^ 



' r;_i m i + m Ti m _i i + 



/! m\ 



(4tt) 



/ + m + 1 (a + c)' +m+1 



0,0,-1 



(4vr) 2 
2a(a + c) 



(74) 



(75) 



A final detail, related to the non-local part of the AA interaction, concerns to the 
implementation of V12 when acting on functions of ri,r 2 and r\ 2 . The operator V12 has 
to be understood as the gradient on the direction r\ 2 = ri — r 2 when the coordinates of 
the center of mass of the AA pair, R cm = (ri + r 2 )/2, are kept fixed. Thus, 

Vu/tri.rj.ru) = j (Vj - V 2 ) /(n,r 2 ,r 12 ), (76) 

= ;(-/--#■ + 2^)/(n,r 2 ,r I2 ), (77) 
2 ari r 2 dr 2 r i2 ari 2 / 

— * 

with Vi(2) the gradient on the direction ri( 2 ) when r 2 (i) is kept fixed. 
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